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Random walks in dynamic random environments and 
ancestry under local population regulation* 
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Abstract 

We consider random walks in dynamic random environments, with an environment 
generated by the time-reversal of a Markov process from the oriented percolation 
universality class. If the influence of the random medium on the walk is small in 
space-time regions where the medium is typical, we obtain a law of large numbers 
and an averaged central limit theorem for the walk via a regeneration construction 
under suitable coarse-graining. 

Such random walks occur naturally as spatial embeddings of ancestral lineages in 
spatial population models with local regulation. We verify that our assumptions hold 
for logistic branching random walks when the population density is sufficiently high. 
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1 Introduction 


Let r]n{x) be a random number of particles located at position x € at time n, where 
■q := {qn)n& '■= iVnix) : x € is a Stationary (discrete time) Markovian particle sys¬ 

tem whose evolution can be described by 'local rules'. We assume that q is in its unique 
non-trivial ergodic equilibrium. Prototypical examples are the super-critical discrete¬ 
time contact process, see l |2.3| l below, or systems of logistic branching random walks, 
see ( |4.4| l in Section 4.1 We consider a random walk X — {Xk)k=o,i,... that moves 'back¬ 
wards' through the medium generated by q, i.e. given 77 , A is a Markov chain and given 
Xk = X, the law of the next increment is a function of 77 in a finite window around the 
space-time point {x, —k). 

Our main result, see Theorem |3.1| in Section provides a law of large numbers 
(LLN) and an averaged central limit theorem (CLT) for X. Very broadly speaking we 
require that the law of an increment of X is close to a fixed symmetric finite-range 
random walk kernel whenever the walk is in a 'good region' and that such good regions 
are sufficiently frequent in a typical realisation of 77 . In particular we assume that 
on suitably coarse-grained space-time scales, the occurrence of good regions can be 
compared to super-critical oriented percolation. The explicit assumptions are rather 
technical and we refer to Sections 13.1113. 2l for details. 

The reversal of the natural time directions of X and q results from and is consistent 
with the interpretation of Xk as the position of the fc-th ancestor of a particle picked 
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from position Xq at time 0. In fact, the spatial embeddings of ancestral lineages in mod¬ 
els with fluctuating population sizes and local regulation are fairly complicated random 
walks in space-time dependent random environments given by the time reversals of 
the local population size processes. In biological applications, they are often replaced 
by ordinary random walks without random environments via an ad-hoc assumption of 
constant local population size; see for example the discussion and references in |[3l and 
Section 6.4 in IfTsTl . 

We verify that in a prototypical discrete spatial population model with local reg¬ 
ulation, namely logistic branching random walks with Poisson offspring distributions, 
the assumptions of Theorem |3.1| are satisfied if the population density in equilibrium 
is sufficiently large. This allows to formulate in Theorem |4.3| a LLN and a CLT for the 
ancestral lineage of an individual sampled from such an equilibrium. Thus, we provide 
at least a partial justification for the aforementioned ad-hoc assumptions from biology 
in the sense that here, an ancestral lineage will indeed behave like a random walk 
when viewed over large space-time scales. This partly answers the question posed in 
lini Chapter 4] in the affirmative. 

As often in the study of random walks in random environments, the main technical 
tool behind our results is a regeneration construction. The details are somewhat in¬ 
volved; in principle, the medium rj can have arbitrary dependence range and in general 
its time reversal cannot be explicitly constructed using local rules. A similar problem 
was faced in lO in the study of a directed random walk on the backbone of an oriented 
percolation cluster. There, the particular structure of oriented percolation allowed to 
jointly construct the medium and the walk under the annealed law using suitable space- 
time local operations (cf. lEJ Sect. 2.1]) and therefrom deduce the regeneration struc¬ 
ture. This construction was extended in II24II to random walks on weighted oriented 
percolation clusters with weights satisfying certain mixing conditions. 

Here, we must use a different approach. Again very broadly speaking, regeneration 
occurs after T steps if the medium r]-T in a large window around Xt is 'good' and 
also the 'local driving randomness' of 77 in a (large) neighbourhood of the space-time 
path {{X^, —m) : 0 < TO < T} has 'good' properties. This essentially enforces that the 
information about r] that the random walk path has explored so far is a function of that 
local driving randomness. Such a time allows to decouple the past and the future of 
X conditional on the position Xt and tj-t in a finite window around that position. A 
difficulty arises from the fact that if a regeneration fails at a given time k, then we have 
potentially gained a lot of undesirable information about the behaviour of r]n at times 
n < —k which might render successful regeneration at a later time £ > k much less 
likely. We address this problem by covering the path and the medium around it by a 
carefully chosen sequence of eventually nested cones, see Figure We finally express 
X as an additive functional of a Markov chain which keeps track of the increments 
between regeneration times and local configurations of rj at the regeneration points. 

Note that random walks in dynamic random environments generated by various 
interacting particle systems, in particular also by the contact process in continuous 
time, have received considerable attention recently; see for example ||ll[2]|4lO[25l|27]|. 
A fundamental difference to the present set-up lies in the time directions. Traditionally, 
both the walker and the dynamic environment have the same 'natural' forwards time 
direction whereas here, forwards time for the walk is backwards time for the medium. 
We also refer to the more detailed discussion and references in lUl Remark 1.7]. 

The rest of this manuscript is organised as follows. We first introduce and study in 
Sectionj^a class of random walks which travel through the time-reversal of the discrete 
time contact process, i.e., 77 is literally a super-critical contact process. We note that 
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unlike the set-up in O, here the walk is also allowed to step on zeros of r]. We use 
this simple model to develop and explain our regeneration construction and obtain a 
LLN and an annealed CLT in the 'p close to 1' regime, see Theorem |2.6[ In Section 
we develop abstract conditions for spatial models and random walks in dynamic ran¬ 
dom environments governed by the time-reversals of the spatial models. Under these 
conditions on a coarse-grained space-time grid we implement a regeneration construc¬ 
tion similar to the one from Section [2] and then obtain a LLN and an annealed CLT in 
Theorem |3.1| In Section we introduce logistic branching random walks, the class of 
stochastic spatial population models mentioned above. An ancestral lineage in such a 
model is a particular random walk in a dynamic random environment, see l |4.10| l. We 
show that this class provides a family of examples where the abstract conditions from 
Section!^ can be verified. We believe that there are several further classes of (popula¬ 
tion) models that satisfy the abstract conditions from Section in suitable parameter 
regions. In Sectionj^we list and discuss such models. 

Finally, we note that a natural next step will be to extend our regeneration con¬ 
struction to two random walks on the same realisation of 77 and to then also deduce a 
quenched CLT, analogous to 0. We defer this to future work. 

Acknowledgements: We would like to thank Nina Gantert for many interesting dis¬ 
cussions on this topic and for her constant interest and encouragement during the pre¬ 
paration of this work. We also thank Stein Bethuelsen for carefully reading a preprint 
version of the manuscript and his helpful comments. Finally, we are grateful to an an¬ 
onymous referee for her or his suggestions that made the presentation more complete. 


2 An auxiliary model 

In this section we prove a law of large numbers and an annealed (averaged) central 
limit theorem for a particular type of random walks in dynamic random environments. 
The model is the simplest and the most transparent among the models that we consider 
in this paper. The proofs here contain already the main ideas and difficulties that we 
will face in the following sections. It will also become clear later how the dynamics of 
ancestral lineages in spatial stochastic population models is related to this particular 
random walk. 

2.1 Definition of the model and results 

We define first the model that generates the dynamic random environment of the ran¬ 
dom walk. Let w := {u}{x,n) : {x,n) G Z‘^ x Z} he a family of i.i.d. Bernoulli random 
variables with parameter p > 0. We call a site {x,n) open if uj{x,n) = 1 and closed if 
uj(x, n) — 0. Throughout the paper |j • || denotes sup-norm. For m <n, we say that there 
is an open path from (x,m) to {y,n) if there is a sequence Xm, ■ ■ ■ ,Xn G Z"^ such that 
Xm = X, Xn = y, \\xk — < 1 for k = m + 1,. ■ ■ ,n and ui{xk, k) = 1 for all A: = m,...,n. 

In this case we write (x,m) {y,n), and in the complementary case {x,m) 74 {y,n). 

For sets A,B <Z Z'^ and m < n we write A x {to} —)■“ B x (nj, if there exist x G A 
and y G B so that {x,m) —>■“ (y,n). Here, slightly abusing the notation, we use the 
convention that uj{x, to) = 1 ^( 2 ^) while for k > m the w{x, k) are i.i.d. Bernoulli random 
variables as above. With this convention for A c Z‘^, m G Z we define the discrete time 
contact process := {r]n)n=m,m+i,... driven by w as 

Vrn and Tj^ (x) . If {Ax {m}—(a;,n)} 5 ^ ^ (2.1) 
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Alternatively = {ri^)n=Tn,m+i,... can be viewed as a Markov chain with and 

the following local dynamics: 


Vn+lix) = 


1 if uj{x, n + 1) = 1 and r]^{y) = 1 for some y s with ||x — y\\ < 1, 


0 otherwise. 


( 2 . 2 ) 


For a distribution y on {0,we write = (y{i)n=m,m+i,... for the discrete time con¬ 
tact process with initial (random) configuration distributed according to p. 

The contact process y^ = (yn)n=m,m+i,... is closely related to oriented percolation. 
In this context A is the set of 'wet' sites at time m and {x e Z'^ : y^ix) = 1} x {n} is the 
n-th time-slice of the cluster of wet sites. Obviously for any p < 1 the Dirac measure 
on the configuration 0 e { 0 , 1 }^ is a trivial invariant distribution of the discrete time 
contact process. It is well known that there is a critical percolation probability pc G (0,1) 
such that for p > pc and any non-empty A c the process y^ survives with positive 
probability. Furthermore, in this case there is a unigue non-trivial extremal invariant 
measure i', referred to as the upper invariant measure, such that, starting at any time 
TO G Z the distribution of y^ converges to as n —oo. 

We assume p > Pc throughout this section. Given a configuration w e {0,1}^ we 
define the stationary discrete time contact process driven by w as 

y — {.yn^nCil, — {^n(^) . X G Z ^ Tl ^ Z} with Pn(^) — ^ X { —oo}—f ^ (a:,n)} ■ (2.3) 

The event on the right hand side should be understood as nm<„{Z‘^ x {to} — {x,n)}. 
In the above notation we have y = y^"^ = Furthermore, since ry is a stationary 

Markov process, by abstract arguments its time reversal is also a stationary Markov 
process with the same invariant distribution; cf. Remark [2.7[ Unless stated otherwise, 
throughout the paper y denotes the stationary discrete time contact process. Many 
other versions of contact processes that will be needed in proofs will be labelled by 
some superscripts similarly to the definition in ( | 2 . 1 | l. 

To define a random walk in the random environment generated by y, or more pre¬ 
cisely by its time-reversal, let 

K := {K„(a;,p) : n G Z, x,y G Z'^} (2.4) 

be a family of random transition kernels defined on the same probability space as y, in 
particular Knix, •) > 0 and Knix,y) = 1 holds for all n e Z and x G Z^. Given 

K, we consider a Z'^-valued random walk X := {Xn)n=o,i.... with Xg = 0 and transition 
probabilities given by 

V{Xn+i = y\Xn = x,k) = Kn{x,y), (2.5) 

that is, the random walk at time n takes a step according to the kernel K„(a;, •) if a: is 
its position at time n. We make the following four assumptions on the distribution of k-, 
see Remark [2. 5 1 for an interpretation. 

Assumption 2.1 (Locality). The transition kernels in the family k depend locally on the 
time-reversal of y, that is for some fixed i?ioc G N 

Kn{x,-) depends only on {w(p,-n), 77_„(y) : ||x - y|| < i?ioc}- ( 2 . 6 ) 

Assumption 2.2 (Closeness to a symmetric reference measure on ry_„(x) = 1). There 
is a deterministic symmetric probability measure Kref on Z'^ with finite range i?ief G N, 
that is Kref(a:) = 0 if ||a:|| > i?ref, and a suitably small Sref > 0 such that 


||K„(a;,a: + •) - Kref( • )IItv < Cref whenever ry_„(a;) = 1. (2.7) 


Here || • ||tv denotes the total variation norm. 
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Assumption 2.3 (Space-time shift invariance and spatial point reflection invariance). 
The kernels in the family k are shift-invariant on x Z, that is. using notation 

■) =i0{z+ ■,m+ ■), 


we have 


Knix, y){uj) = Kn+m{x + Z,y + z){d''’"^Uj). (2.8) 

Moreover, if g is the spatial point reflection operator acting on w, i.e., Quj{x,n) = 
a;{—X, n) for any n G Z and x G Z'^, then 

/c„(0,2/)(w) = K„(0,-j/)(pw). (2.9) 

Assumption 2.4 (Finite range). There is € N such that a.s. 

Knix,y)=0 whenever \\y — x\\ > Rk- (2.10) 

Remark 2.5 (Interpretation of the assumptions). The Assumptions |2 . 1}f2 .4| are natural 
as we want to interpret the random walk as the spatial embedding of an ancestral lin¬ 
eage in a spatial population model, in which roughly speaking, children (if any present 
at a site) choose their parents at random from a finite neighbourhood in the previous 
generation. See also Section]^ and in particular the discussion around ( |4.10| l. 

By ( |2.5| l and ( |2.6| l, we can, and often shall think of creating the walk from rj and oj in 
a local window around the current position and additional auxiliary randomness. 

The main result of this section is the following theorem. Its proof is given in Sec¬ 
tion [231 

Theorem 2.6 (LLN and annealed CLT). One can choose 0 < Cref sufficiently small and 
p sufficiently close to 1, so that if k satisfies Assumptions |2.l|[2.4| then X satisfies the 
strong law of large numbers with speed 0 and an annealed (i.e. when averaging over 
both w and the walk) central limit theorem with non-trivial covariance matrix. A corres¬ 
ponding functional central limit theorem holds as well. 

Remark 2.7 (Time-reversal of ij, oriented percolation interpretation). In ||5j|, the station¬ 
ary process 77 was equivalently parametrised via its time reversal 

^ := {^(x, n) : X G Z*^, n G Z} with ^(x, n) := 77_„(x). 

Then ^ is the indicator of the backbone of the oriented percolation cluster and it was 
notationally and conceptually convenient to use in O, not least because then the me¬ 
dium ^ and the walk X had the same positive time direction. 

Here, we keep ry as our basic datum because we wish to emphasise and in fact later 
use in Section the interplay between the medium rj, interpreted as describing the 
dynamics of a population, and the walk X, cf. I |2.5| l above, describing the embedding 
of an ancestral lineage. Furthermore, in the more general population models, as the 
one studied in Section [^ for instance, there will be no natural parametrization of the 
time-reversal of 77. 

Note that the assertions of Theorem |2.6| are in a sense conceptual rather than prac¬ 
tical because the proofs of preliminary results in Section |2.3| require 1 — p to be very 
small. Situations with p > pc but also I — p appreciably large require an additional 
coarse-graining step so that the arguments from Sectioncan be applied. 

In order to prove Theorem |2.6| we will construct suitable regeneration times and 
show that the increments of these regeneration times as well as the corresponding 
spatial increments of the walk have finite moments of order b for some b > 2. This 
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construction is rather intricate. The main source of complications stems from the fact 
that in order to construct the random walk X one should know uj and r? in the vicinity of 
its trajectory; cf. Remark [2.51 above. While it is easy to deal with the knowledge of w's, 
because they are i.i.d., the knowledge of 77 leads to problems. Due to definition ( |2.3| l of 
77 and Assumption |2.1| on k, this knowledge provides non-trivial information about the 
past behaviour of 77 and therefore also about the future behaviour of X. Both is not 
desirable at regeneration times. 

More precisely we need to deal with two types of information on 77 . The first type, 
the negative information, that is knowing that rjnix) = 0 for some n and x is dealt 
with similarly as in lO. The key observation is that such information is essentially 
local: To discover that Pni^x) = 0 one should check that x {— 00 } {x,n) which 

requires observing uj's in a layer x {n — T,... ,n} where T is a random variable with 
exponentially decaying tails. The second type of information, the positive one, that 
is knowing that 77 „(x) = 1 , is removed by making use of strong coupling properties 
of the forwards-in-time dynamics of 77 . When at time —t we have rj-t > 
pointwise, then there is a substantial chance that every infection, i.e., every of 77 , 
inside a growing space-forwards-time cone with base point (x, —t) can be traced back 
to (x + X {—t}. Furthermore, whether this event has occurred can be 

checked by observing the restriction of p-t to x + {—T,... ,LY and the w’s inside a 
suitably fattened shell of the cone in question, in particular without looking at any Pmiy) 
for m < —t; see ( |2.27| ) and Lemma [2 .13 1 below. We will construct suitable random times 
T at which this event occurs for 77 at the current space-time position {Xt, —T) of the 
walker and in addition the space-time path of the walk up to T, {{Xk^—k) : 0 < k < T}, 
is completely covered by a suitable cone. Such a time T allows to regenerate. 

For the proof of Theorem |2.6| we first collect some results on the high density dis¬ 
crete time contact process in Section |2.2[ We then rigorously implement the regen¬ 
eration construction sketched above in Section |2.3| Finally we prove Theorem |2.6| in 
Section [2~4l 

2.2 Some results about the contact process 

This section contains several estimates for the discrete-time contact process 77 that will 
be crucial for the regeneration construction. The main results of this section are the 
estimates given in Lemma [2.11| and Lemma [27^ We start by recalling two well known 
results. 

Lemma 2.8. Forp > pc there exist C{p),c{p) G (0,(x;) such that 

P(Z‘^ X {-n} (0,0) and x {- 00 } 74 “ (0,0)) < C'(p)e-‘=(P)”, uGN. (2.11) 

Moreover, we have limsuppyi^ C{p) < 00 and limpy<i c{p) = 00 . 

Proof. Due to self-duality of the contact process this is a reformulation of the fact that 
for p> Pc and 77 ^°^ = there exist C'(p), c(p) € ( 0 , 00 ) such that 

P( 77 );°^ ^ 0 and 77 ^°^ eventually dies out) < C'(p)e“'^^^^”, tt G N. (2.12) 

For a proof of the latter assertion we refer to e.g. IIIOII18II or Lemma A.l in |[5l. □ 

Lemma 2.9. Let 77 ^°^ = (? 7 n°^)n=o,i.... and let p'' = {Pn)n=o.i,...- where v is the upper 
invariant measure. For p sufficiently close to 1 there exist Scoupi > 0, C < co, c > 0 such 
that 


P(77f >(x) = 77(((x) for all ||x|| < Scoupin | pi^^ ^ O) > 1 - neN. (2.13) 
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Figure 1: Possibly overlapping finite clusters starting at F = {(xi,ti)} that appear in 
Lemma [2. 11 [ Here k = 4, Di — 2, D 2 = ^,03 = 3, D 4 — 1, hence M = 3, Si = 1, S 2 = 
3 ,53 = 6. 

Proof. For the contact process in continuous time, this is proved in 01211 . see in par¬ 
ticular (33) and (34) in Proposition 6 there. Although literally, IfT^ Eg. (34)] refers to 
conditioning on { 77 ^°^ survives} the result follows in view of \2.12\ . □ 

Remark 2.10. In 0160 it is shown (literally, for the contact process in continuous time) 
that for any p > Pc and a > 0, there is a (7 < cx) such that the probability on the left hand 
side of ( |2.13| l is bounded below by 1 — Cn~°-. 

More recently, in 0170 large deviations for the continuous time contact process in a 
random environment were studied. Among other results it is shown that exponential 
decay as in l |2.13| l holds in the supercritical case; see Theorem 1, Eg. (2) there. 

The first main result of this section is the following lemma on controlling the prob¬ 
abilities of certain negative events; cf. Lemma 7 in Section 3 of QUO for a related result. 
Lemma 2.11. For p large enough there exists s{p) G (0,1] satisfying liinp^ie{p) = 0 
such that for any V = {{xi,ti) : 1 < i < k} c x Z with ti > t 2 > ■ ■ ■ > tu, we have 

P( 77 t(a;) = 0 for all (a;, t) gV) < e{p)^. (2.14) 

Remark 2.12. In our proof of l |2.14| l it is essential that all tfs are distinct. For a general 
set V c X Z space-time boundary effects can play a role so that the decay will only 
be stretched exponential in \V\. For a concrete example in the case d = 1 consider the 
set 


V = {{x,n) gZ X Z : |a;| < h,0 < n < h} 

and the event 

B = {oj{y, 0) = 0, |j/| < /i -I- 1} n {uj{±{h + 1), /c) = 0, fc = 0, 1,..., h}. 

Then, obviously we have B c {r]n{x) = Ofor all {x, n) G V} and P(i?) = (1 — 

Note however that in Corollary 4.1. in 02311 it is shown that the upper invariant meas¬ 
ure of the continuous time contact process on Z'^ dominates stochastically a product 
measure on {0, 1}^ . Therefore a bound analogous to ( |2.14| l does hold in the situation 
when all t/s are egual. Lemma [2.11| can be seen as a space-time extension of that result 
in the discrete time case and in the 'p large enough’ regime. 

Proof of Lemma [2.11i An immediate conseguence of Lemma [T8| is that for every p > p^ 

P(Z’'x (0,0) and Z’'x {-c3o} (0,0)) < n = 0,l,2 ,..., 

(2.15) 
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with some ci = ci{p) > 0 satisfying limp/.ici(p) = oo. To prove l |2.14| l, we extend the 
finite sequence {ti,... ,tk} to an infinite sequence via tk+j ■= tk— j, j = 1,2 ,..., and put 

A := min {£ G N : x {U+i} [xi^U)] . (2.16) 

Note that the random variables Di are upper bounds on the heights of the backwards- 
clusters of open sites attached to {xi,ti) given by (see Figure 0 

{(y, to) G Z'' X Z : TO < U, {y, m) (xi, U)). 

For each {xi,ti) G V we have r]t-{xi) = 0 if and only it Di < oo. Thus the left-hand side 
of ( |2.14| l satisfies 

k 

F{r]t{x) = 0 for all {x,t) G V) =p(^P|{A < oo}). (2.17) 

On the event < oo} we further define 

S*! = 1,5'2 = 5*1 + Ds ^) • • ■) >S'i+i = Si + Dg. as long as Si < k, 

and let M be such that S'm-i < k < Sm: see Figure[^ For z = 1,..., M we set A := Dg^ 
and Di := oo for i > M. Finally we set 

I{m, k) = {{di ,..., dm) G N"* : di -I- • • • -|- dm-i "£ k < di + ■ ■ ■ + dm}, (2.18) 

and for (di,..., dm) G l{m, k) we write 

m(1 ) = 1, u(2) = u(l) + di, ..., u(to) = u{m — 1) -I- dm-i- 

Then we have 

k k 

f|{A<00}cU y {5i =di,...,Dm = dm}. (2.19) 

i—1 'm—l(^di,...^dm)^^{Di,k) 

Note that 

— C?x, . . . , Dm — ^ 

m 

= y ^ {tu{j)+dj-l} i^uU)du{j))} O X {tu{j)+dj} i^u{j)du{j))} ■ (2.20) 

i=i 

The events in the intersection on the right-hand side depend on co restricted to disjoint 
sets and are thus independent. Furthermore we observe that the event 

|Z X ^ (.^u{j)jtu(j))} O ^ {^u{j)+dj} i^U'{j):tu{j))} 

enforces that {xu(j),tu(j)) is the starting point of a finite (backwards) cluster of height 
at least — t„(j) > d^ — 1 (when dj = 1 this means that uj{xu(j),tu{j)) is closed, 

which also gives a factor 1 — p < 1). Hence, using l |2.15| l we obtain 

P(i5l = di, . . ..Dm = dm) 

m 

= JIP^Z^* X {tu(j)+d,-l} {Xu{j),tuij)),'^'‘ X {tu(j)+d,} 7^“^ iXu{j),tuU))}) 

i=i 

m 

< _ g-ci(p) EJLi dj ^ 

i=i 

(2.21) 
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Now ( |2.21| l with ( |2.19| l imply 


F{r]t{x) = 0 for all {x,t) &V) < E E 


=-ci(p) ETLi dj 


m—1 (di. .,dJ 7 ^)GX(m,fc) 


= ^ ^ e-^i(P>-#{(di,...,d„)eI(m,fc):di + ---+d„=s}. (2.22) 

m—1 s—k-\-l 

By definition of I{m, k) for s > k + 1 we have 

; ^m) € ^{rn^ k'j . di -t- * * * -t- dm — s} 

= ■ j dm) G dl{rn^ k) : di dm = fc + 1} = 


/c 

m — 1 


Thus, the right hand side of ( |2.22| l can be bounded by 

^ p-ci(p)(fc+l) 

2'^y y = k2^^ - 

^ ^ 1 — e“^i(p) 

m—1 s—k-^1 


<2^. (2.23) 


(2.24) 


yielding the claim of the lemma. 


□ 


The second main result of this section. Lemma [2.131 below, is the crucial tool in the 
construction of a certain coupling which will be useful to forget the positive information 
about r] in the regeneration construction. To state this lemma we need to introduce 
more notation. For A c let rj^ = (? 7 ;^)„=o,i,... be the discrete time contact process 
as defined in iD- For positive b, s, h we write (denoting by Z+ the set non-negative 
integers and by 11-112 the £^-norm) 

cone(5, s,/i) := {(x, n) S X Z+: ||x ||2 < & + sn, 0 < n </i}. (2.25) 

for a (truncated upside-down) cone with base radius b, slope s, height h and base point 
(0,0). Furthermore for 


(^inn ^ ^out UUd Sinn ^ Souti (2.26) 

we define the conical shell with inner base radius 6 inn, inner slope Sinn, outer base 
radius 6 out, outer slope Sout, and height /iGNU{oo}by 

Os(&inn: ^out; Sjnn; Sout: — { (^5 Xl) G Z X Z . ^inn H“ ^inn^ < ll^^lb < ^'out + SoutUjO <n<h]. 

(2.27) 

The conical shell can be thought of as a difference of the outer cone cone(bout,Sout,h) 
and the inner cone cone( 6 inn, Sinn, h) with all boundaries except the bottom boundary of 
that difference included; see Figure 

Let 77 “ := (? 7 n )n=o,i.... be the contact process in cs( 6 inn, (>out, Sinn, Sout, oo) with initial 
condition 77 §=(x) = l{&.„„<|| 2 :|L<b„„t} and 

{ 1 if (x, n -f 1 ) G Cs( 6 inm ^out; Sjnn; Sout; Oo), U;(x, U -f 1 ) — 1 
and ? 7 “( 2 /) = 1 for some y S Z"^ with ||x - y\\ < 1 , 

0 otherwise. 

We think of 77 “ as a version of the contact process where all w's outside the conical 
shell cs(6inn, (>out, Sinn, Sout, 00 ) have been set to 0. We say that 77 '^® survives (in all parts 
of the conical shell) if for all tt S Z+ there is x G Z'^ with = 1. In the case 


9 













Figure 2 : The left figure shows cone( 5 , s, h). The grey region in the figure on the right 
(without the bottom line) shows cs( 5 inn, &out, Smn, Sout,h) in x Z. 

d = 1 we require additionally that for all n S Z+ there \s x & Z+ and y S Z_ with 
= dniy) = 1 - (Here Z_ denotes the non-positive integers.) For a directed path 

7 = ((a;m,m), (xm+i,TO + 1),..., (a;„,n)), m < n, 0;^ S with - XiH < 1 

( 2 . 28 ) 

we say that 7 crosses the conical shell cs(6inn, ('out, Smn, Sout, 00) from the outside to the 
inside if the following three conditions are fulfilled: 

(i) the starting point lies outside the outer cone, i.e., ||a;m||2 > ('out + msout, 

(ii) the terminal point lies inside the inner cone, i.e., ||a:„||2 < ('inn + nsinn, 

(hi) all remaining points lie inside the shell, i.e., G cs(6inn,('out, Sinn, Sout,00) for 

i = m + 1 ,... ,n — 1 . 

We say that 7 intersects 77 '^^ if there exists i G {m -I-1,..., n — 1 } with r]f{xi) = 1 . Finally 
we say that 7 is open in cs(5inn, ^out, Smn, Sout, 00) if oj(xi, i) = 1 for alH = to + 1,..., n — 1. 

Note that if in the case d = 1 the process 77“ survives in a ray of a conical shell and 
7 is a path that is open in this ray then by geometric properties of directed paths and 
clusters 7 necessarily intersects 77“. In the case d > 1 however, even if 77^^ survives, it is 
in principle possible that an open path crosses the conical shell cs(()inn, ('out, Sinn, Sout, 00) 
without intersecting 77“. The next lemma states that the probability of that can be made 
arbitrarily small. 

Lemma 2.13. Assume that the relations in ( | 2 . 26 | l hold and consider the events 
Gi := {77'^^ survives}, 

G2 := {every open path 7 that crosses cs(6inn, ('out, Sinn, Sout, c») intersects 77^^}. 

For any e > 0 and 0 < Sinn < Sout < 1 one can choose p sufficiently close to 1 and 
('inn < ('out Sufficiently large so that 

P(GinG2) > 1-e. 

Remark 2.14 (Observations concerning Gi and G2). The meaning of the event Gi is 
clear. Let us just note that it is essential that the relations in l | 2 . 26 | l hold. In particular, 
in the case Sjnn = Sout survival of is only possible in the trivial case p = 1 . 

To understand the importance of the event G2, observe that if a path 7 as in ( | 2 . 28 | l 
crosses and is open in cs(()inn, ('out, Sinn, Sout, 00), and also intersects 77'=^ then necessarily 
Vnixn) = n) for the terminal point (x„, n) of the path. Thus, on Gi (7 G2 the values 
of the contact process inside the inner cone, that is for {x,n) with ||x||2 < (>inn + nsinn, 
are independent of what happens outside of the shell; cf. I | 2 . 59 | l and Lemma [ 2.221 
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Proof of Lemma f2.13l The proof consists of two steps. In the first step we prove the 
assertion for the case d = 1. The second step then uses the assertion for d = 1 to give a 
proof for d> 2. 

Throughout the proof of this lemma for r > 0 and x G we denote by Br{x) the 
closed ball of radius r around x, i.e., Br{x) = {y : \\x — y \\2 < r}. 

Step 1. Consider the case d = 1. We first check that the discrete time contact process 
survives with high probability in any oblique cone when p is large enough. To this end 
for 0 < Si < S2 < 1 and 6 e N we set 


C&,si,s 2 ■= {(x,n) : X G Z, n G Z_|_, sin < x < s^n + &}. 


Furthermore we let f] := (77„)„=o,i,... be the discrete time contact process in Ct,_si,s 2 
starting from % = ]l[o,6]nz and with all w's outside set to 0. Finally set Xn = 

\b/2 + n(si + S2)/2J, r„ = n(s 2 - si)/4. 

Claim 2.15. (i) For every 0 < si < S 2 < 1 and £ > 0 there is h large and po < 1 such that 
for all p> po, fj survives with probability at least 1 — e. 

(ii) Moreover, there exist c,C G (0,oo) so that on the event {iy survives}, with prob¬ 
ability at least 1 — restricted to the ball Br„{xn) can be coupled with the 

unrestricted process = {r]n)n=o,i,... started from the upper invariant measure i/, that 
is fjn{x) = T]n{x) for all X G Br^{Xn). 


The claim (i) follows using the same arguments as in the proof of Theorems 1 and 2 
of 0 , where an analogous statement is proved for continuous time contact process in a 
wedge. Moreover, 0 use for the proofs a coarse-graining construction and comparison 
with oriented percolation. That links the p roble m in the wedge with a suitable shift 
of the contact process 77 ^°^ used in Lemma 2.9 Combining that lemma with coarse- 
graining then yields the claim (ii). 

As noted in the paragraph above Lemma 2.13 we have Gi c G 2 in the case d = 1. 
Thus, by using the above argument for the oblique cone twice we see that the assertion 
of the lemma holds in the case d = 1 . 


Step 2. Consider now d > 1. Since the probability of Gi is increasing in dimension, it 
can be bounded from below by the same reasoning as in d = 1. It remains to show that 
the probability of can be made small by choosing 6 inn, &out and p appropriately. 

To this end for n > 0 we set d(n) = (6inn + &out)/2 + n(sinn + Sout)/2 and define 

Mn ■= {x gZ‘^ : ||a ;||2 G [d(n),d(n) + 1]}. 


Furthermore for x G Mn we consider the event 


Bn{x) ;= {rjnix) = 0 and {x,n) contained in an open path 7 crossing the conical shell}. 
Finally, we define the backward cluster of {x, n) by 

bC(x,n) ;= {(y,m) G X Z : m < n, {y,m) —(a;,n)}, 


and for m <n we set 


hCm{x,n) ■■= {y gZ‘^ : {y,m) G bC(a;,n)}. 


Assume that Bn{x) occurs. Then, by trivial geometrical arguments there is a small 
constant p depending on the parameters of the shell so that for M = [(1 — p)n\ the set 
hCM{x,n) is not empty. Moreover, self-duality of the contact process and Lemma 2.9 
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imply that there exist Scpi > 0 (here we think of Scpi ~ ScoupiP for Scoupi from Lemma [2^ 
and c ^ 0 such that 13 x cz cs(6inn; (^out? "^inn; ^out 5 oo) and 


with probability bounded below by 1 — e“™, the indicator function of 

the set bCM(a;,n) can be coupled inside Bs^^^nix) with the set of I’s (2.29) 

under the upper invariant measure v of the (full) contact process. 

Fix p large enough so that the density of I's under v is strictly larger than 1/2; this 
is not a restriction in the parameter region that we consider. On the one hand this 
requirement means heuristically that the set bCM(a;, n) n {Bs^^^n{x) x {M}) is large with 
high probability. Thus we must have — 0 for oil {z,M) in this set. On the other 

hand, using d = 1 -arguments we will show that this is not possible. 

To this end, depending on the previous parameters, we fix 5 > 0 small, and unit 
vectors u® G 1 < i < N, with N sufficiently large so that for every x G Mn, there is 
ani < N such that 


the length of the intersection of the half-line {tu* : f > 0 } with the 3 qj 

(real) ball {y : ||a; — y \\2 < Scoupi^r} has length at least 6 n. 

Observe that N and u*'s can be chosen independently of n. For i < N, let a* = 

be a self-avoiding nearest neighbour path in Z‘^ approximating the half-line {tv^ : t> 0 } 

given by 


(a) Off) = 0, 

(b) Of* makes steps only in direction of u*, that is for every coordinate k = 1 ,... ,d and 

j > 0 one has + 1) - > 0, 

(c) a* stays close to fn®, that is {a® : j G N} C {to® + z : t > 0, z G ||z|| < 2}. 

Using ( |2.29| l, ( |2.30| l and large deviation estimates for the density of I's under 
see Thm. 1] (literally, proved there for the continuous-time contact process) we see 
that there is c > 0 such that that for i < N satisfying l |2.30| l, 

W n n bCM(x,n)| 

I a® n Bs^^,n{x)\ 



Bn{x) ) > 1 — ce 


(2.31) 


We now use the result of Step 1 and the last claim to bound the probability of Bn{x). 
To this end we define contact process 77 ^®) as the contact process restricted to the set 
W* ;= (a® X Z+) n cs(&inn, ('out, Sinn, 00) Started from lw’n(Z‘*x{o})- Observe that W® con¬ 
tains an isomorphic image of Cb^sus 2 > for some 6 , si,S 2 - Thus, the contact process ryf®) 
dominates a corresponding image of a contact process in Cb^si,s 2 - 

Let S defined be the event {Ty^®^ survives for every i < N}. Given e > 0, choose e' so 
that (1 — e')^ > 1 — s/2. Using the result of Step 1 with e' replacing e, we see that 


P(5)>l-e/2. (2.32) 

Moreover, on S, for suitable x^/^ with probability at least 1 — , Ty^®^ can be coupled 

io Brnf ix\J) X {M} to a stationary contact process. Moreover, the parameters can be 
chosen so that for every x G Mn, there is x^/^ for which Bs^^^„{x) n a® C Br^{x^/^). 
Hence, using again large deviation arguments for the density of I's we obtain 


\{y G a® n : ?y^^( 1 /) = 1}| 


la® n B, 


*(2;) I 


> 1/2 


5 > 1 — ce 


(2.33) 
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Comparing ( |2.31| l with the last display, we see that for every x € Mn, 


P(e„(a;)|5) < e-^". (2.34) 

Assume that occurs. Since there are at most polynomially many x € Mn, there must 
be {x,n) e x {£}) so that Bn{x) occurs. It follows that for p sufficiently large 

we have 

P(G2|5) > l-e/2 (2.35) 

and therefore P(Gi n G 2 ) > 1 — e, as required. □ 

2.3 Regeneration construction 

In Theorem |2.6| we claim that the speed of the random walk X is 0. As an intermediate 
result we will show that the speed is bounded by a small constant. This will be needed 
for the regeneration construction. 

Lemma 2.16 (A priori bound on the speed of the random waik). If the the family of 
kernels k satisfies Assumptions |2.1H2.4| then there are positive finite constants Smax. c 
and C so that 


P(||X„|| > wn) < Ge-™, n e N, (2.36) 

in particular limsup„_>oo||X„||/n < Smax almost surely. The bound Smax can be chosen 
arbitrarily small by taking Sref "C 1 (where Cref is from Assumption [2]^ and 1 — p <C 1. 

Proof. With the percolation interpretation in mind, we say that a space-time site {x, k) 
is wet if 'qk{x) = 1, and dry if r]k{x) = 0. Let r„ be the set of all n-step paths 7 on 
starting from 70 = 0 with the restriction ||7i — 7^-1 1 | < R ^, i = I ,... ,n, where is the 
range of the kernels from Assumption |2.4[ For 7 G r„ and 0 < < (2 • • • < ffc < n we 

define 


'■= id-eile) = 0 for alH G {ii,..., 4}}, 

■= iv-d'li) = 1 for alH G {1,..., n} \ {fi,..., ifc}}- 


Let Hn ■= #{0 < i < n : r]-i{Xi) = 0} be the number of dry sites the walker visits up to 
time n and set K := max^gzd{«:ref(x)} + Cref- For fc G {1,..., n} by Lemma [2. 11 we have 


V{Hn = k)= Y. E 

^ E E 

0<ii<‘“<ik<n 7 GrTi 


(2.37) 


It follows 


n 

V{Hn >Sn)< Y 

k—[nS\ 




k 


00 

<{2RiKY Y 

k— [nE 


{2RiKY 


{e{p)/KY 

1 - s{p)/K 




(2.38) 


with Cl, C 2 G (0, 00 ), when 5 > 0 is sufficiently small and p> po = PoiS, Cref)- 


13 










Writing = (X„ i, ..., we can couple the first coordinate (-Y„,i)n=o,i,... of the 
random walk X with a one-dimensional random walk X — {Xn)n=o,i,... with transition 
probabilities given by 

^n— 1 — — (1 ^ref) ^ ^ ^ref (O? ^2^ • ■ • “t“ £^ref (^) ; X € ^ 

(x 2 ,...,xd)eZ'^~^ 

(i.e., X takes with probability 1 — £ref a step according to the projection of K^ei on the 
first coordinate and with probability e^ef simply a step of size to the right) such that 
for all n e N 


^ ^n — Hn “f Rn^n- 


Then, we have 


Ln<5J 

P {Xn^i > sn) < P {Hn > 5n) + P (X^ i > sn, = k) 

fc=0 

[nij ^2 

< P {Hn > 5n) + ^ P yXn-k > sn - 

k=0 

< P {Hn > 6n) + (5nP > (s — dRK)n'^ . 

The estimates ( |2.38| l, ( |2.39| l and standard large deviations bounds for X show that 


P {XnA > sn) < holds for all n € N (2.40) 

with C 3 , C 4 S (0, oo) when s — SR^ > E[Xi] = £reff?«;- By symmetry we have an analogous 
bound for P {Xn,i < —sn). The same reasoning applies to the coordinates Xn,2, • ■ ■ , Xn,d- 
Thus, we have 

d 

P(||X„|| > sn) < ^P(|X„,,| > sn) < 2dc3e-"^”. (2.41) 

i=l 

In particular we have limsup„_>ooll-’^n||/n < s a.s. by the Borel-Cantelli lemma. □ 

Denote the i?ioc-tube around the first n steps of the path by 

tube„ := {{y, -k) ■■ Q < k < n,\\y - Xk\\ < i?ioc}- (2.42) 


For (x,n) G Z‘^ x Z let £{x,n) be the length of the longest (backwards in time) directed 
open path starting in {x, n) with the convention £{x, n) = — 1 if w(x, n) = 0 and £{x, n) = 
00 if r]n{x) = 1. For each (x, n) we define its determining triangle by 


D{x, n) 


0 , 

{{y,m) : \\y - x|| < (n - m), n - £{x,n) 


1 < m < n}, 


if r]n{x) = 1, 
if 77„(x) = 0. 


(2.43) 


The idea is that if ry„(x) = 0, i.e. (x,n) is not connected to Z‘^ x {—oo}, then this inform¬ 
ation can be deduced by inspecting the w's in D{x,n). By definition of £{x,n), in this 
case there must be a closed contour contained in D{x,n) which separates (x,n) from 
Z‘^ X {—oo}; see Figure]^ Note in particular that D{x, n) = {{x, n)} if uj{x, n) = 0. 

When constructing the walk X for n steps we must inspect w and y in tube„ (cf. Re¬ 
mark [2]^. By the nature of y, this in principle yields information on the configurations 
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{x, n) 



Figure 3: A caricature of the determining triangle D{x,n) with a closed contour. The 
height of the triangle is i{x, n) + 1. 


r]_k, k > n that the walk will find in its future. Positive information of the form rj^iy) = 1 
for certain m and y is at this stage harmless because y has positive correlations and in 
view of Assumption |2.2| this suggests a well-behaved path in the future. On the other 
hand, negative information of the form rj^iy) = 0 for certain m and y is problematic be¬ 
cause this increases the chances to find more O's of 77 in the walk's future. In this case 
As sumption 1 2.2 1 is useless. In order to 'localise' this negative information we 'decorate' 
the tube around the path with the determining triangles for all sites in tube„ (obviously, 
only zeros of ry matter) 

dtube„ = y D{y,k). (2.44) 

(y,/c)Gtube„ 

Define 

Dn ■=n + max {£{y, -n) + 2 : ||X„ - y\\ < R\oc, £{y, -n) < 00 }. (2.45) 

Note that is precisely the time (for the walk) at which the reasons for 77 _„(j/) = 0 for 
all y from the i?ioc-neighbourhood of Ai„ are explored by inspecting all the determining 
triangles with base points in BR^^^{Xn) x {—n}. The information rj-niy) = 0 does not 
affect the law of the random walk after time Note that the 'height' of a non-empty 
triangle D{y, —n) is £{y, —n) + 1. This is why —n) + 2 appears in definition l |2.45| l. 

Now, between time n and Dn the random walk might have explored more negative 
information which in general will be decided after time and will affect the law of the 
random walk thereafter. To deal with this cumulative negative future information we 
define recursively a sequence 

(To := 0, ai '■= min <m> (Ti_i : max < m >, z > 1. (2.46) 

L (7i_i<n<m J 

In words, at is the first time m after m-i when the reasons for all the negative informa¬ 
tion that the random walk explores in the time interval ai-i ,..., m are decided 'locally' 
and thus the law of the random walk after time ai does not depend on that negative 
information. The a^ are stopping times with respect to the filtration T = (J^„)n=o,i, 2 ...., 
where 

Xn ■■= cr{Xj : 0 < j < n) V a{r]j{y),uj{y,j) : {y,j) S tube„) V a{uj{yJ) : {y,j) S dtube„). 

(2.47) 


Note that by construction we have r]-^.{y) = 1 for all y e 

Lemma 2.17. When p is sufficiently close to 1 there exist finite positive constants c and 
C so that 


P((Ti+i - (Ji > n I F„.) < Ce ™ 


for all n = 1,2,..., z = 0,1,... a.s., 


(2.48) 
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in particular, all tJi are a.s. finite. Furthermore, we have 

foralH = 0,l,... a.s., (2.49) 

where V' denotes stochastic domination. 

Proof. Throughout the proof we write := (2i?K + 1)'^ and i?ioc := (2i?ioc + l)'^ for the 
number of elements in respectively in 

Consider first the case f = 0 in ( |2.48| l. The event {tJi > n} enforces that in the 
i?ioc-vicinity of the path there are space-time points (%,—j) with rj-j{yj) = 0 for j = 
0,1,..., n. For a fixed choice of the yj's by Lemma [2. 11| the probability of that event is 
bounded by £(p)”. We use a relatively crude estimate to bound the number of relevant 
vectors (j/o, 2/i, • ■ • > Vn) C as follows. There are i?” possible n-step paths for the 

walk. Assume there are exactly k time points along the path, say 0 < mi < • • • < m^ < n, 
when a point {yrm.—rrii) G x {—m*} with = 0 is encountered and 

hence the corresponding 'determining' triangle D{ymi,—mi) is not empty (when n > 1, 
we necessarily have mi = 0 or mi = 1). 

For consistency of notation we write mk+i = n. Then the height of D{ymi,—mi) is 
bounded below by m^+i — m^. For a fixed n-step path of X and fixed mi < • • • < m^, 
there are at most many choices for the i = 1,..., fc, and inside D{y^.^ —rrii) 
we have at most choices to pick yjni+i,ymi+ 2 , ■ ■ ■ ,ymi+i-i (start with ymt, 

then follow a longest open path which is not connected to x {—oo}, these sites are 
necessarily zeros of y). Thus, there are at most 

n k n 

K E E ^loc n = i?” E 

/l‘= 1 mi<---<mfc<mfc + i=n i—1 k—1 


itioc-ft, 


i—k 


^RZiRioc + R.y 


possible choices of {yo,yi,... ,yn) and hence we have 

P(cri > n) < {Rk{Rioc + Rk)£{p)Y- 


The right hand side decays exponentially when p is close to 1 so that e(p) is small 
enough. For general i > 0 l |2.48| l follows by induction, employing l |2.49| l and the argu¬ 
ment for i = 0. 

In order to verify ( |2.49| l note that the stopping times Oi are special in the sense that 
on the one hand, at a time Oi the 'negative information' in , that is the knowledge of 
some zeros of y in the i?ioc-neighbourhood of the path, has been 'erased' because the 
reasons for that are decided by local information contained in Ren ■ On the other hand, 
the 'positive information', that is the knowledge of ones of y, enforces the existence of 
certain open paths for the uj's. And this information is possibly retained. Thus, ( |2.49| l 
follows from the FKG inequality for the w's. □ 


Corollary 2.18 (Reformulation of Lemma 2.11 . For any F ..., (xkRk)} C 

X N with ti < t2 < ■ ■ ■ < tk and e(p) as in Lemma 


2.11 


we have 


=0forall {x,t) G 14| J",..) < £{pf . (2.50) 

Proof. The assertion is an easy consequence of ( |2.49| l and Lemma [2.11[ □ 

For f G N we define Rt ■= inf{i G Z+ : ai> t} and for m = 1,2,... we put 



else. 


(2.51) 


In words, is the length of the time interval {ui-i,ai] which contains t and Tm is the 
length of the (m — l)-th interval before it. 
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Lemma 2.19. When p is sufficiently close to 1 there exist finite positive constants c and 
C so that for all i, n = 0,1,... 

> n I < Ce“™ a.s. on {tJi < f}, (2.52) 

and generally 

>i + m, > n I for m = 1, 2,... a.s. on {ai < t}. (2.53) 

Proof. For l |2.52| l, we have 

t-1 

= P(cri+i > f V (n + CTi) I J'o-i) + ^ ^ P(crj = £, (Tj+i > f V (£ + n) I J'cr,) 

j>i l^(7i + l 

t-1 

< Ce-™ + ^ j>i : aj=£\F^^) 

t—n—1 t—1 

< ^ + !{„+!<*} ^ < C(l + + n)e-^" 

+ l i—t—n 

where we used Lemma r2 .17 1 and 

P((Tj = f, CTj + l > f V (£ + n) I J'crj = E[l{^^=£}P(tT^ + l - cr^- >(<-.() V n I J'crj 
in the first inequality. 

Similarly for m > 2 (we assume implicitly that cr^ <t — n — m—1 for otherwise the 
conditional probability appearing on the right-hand side of l |2.53| l equals 0) 

t—m—n t—m-\-l 

P(r^^ > n I Fai) = X! X! X! ^ ^ ^ ^ I 

t—m—n t—m-\-l 

j >2 fc=tTi + l £—k-\-n 

t—m—n t—771+1 

< ^(to - 1) XI Yl X 

fc=(7i + l ^=fc+n j^i 

t—m—n t—771+1 

<C'^(m-l) X X exp(-c^£) 

A;—tTi + 1 £—k-\-n 

where we used in the first inequality that 


j+m 

{u'j + X f, < f, (T^-|_772 ^ f} ^ (Jy—l ^ m — l } 

r=j+2 

together with Lemma |2.17| and then argued analogously to the proof of l |2.52| l for the 
second inequality. For m = 2, the chain of inequalities above yields the bound 

t — n oo 

P(f+ > n I X (i- 

k=CTi + l 1=0 
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whereas for m> 2 we obtain 


r(t) 


> n I < C'^{m — 1) ^ 


. ,,, exp ( — + n)) 

^ck — ctj (m —1) 


k—ai + 1 


C'^{m — 1) 


t — m—n 


< 


C‘^{m — 1 ) 


k—{Ti-\-l 
^c{t-n)/{m-l) 


ET exp ( - c^n - c;;^) 


C'^{m — 1) e“™ 

1 _ - 1 


gc/(m-l) _ 

< C"(to- l)2e-™. 


Thus, ( |2.53| l holds (with suitable adaptation of the value of the prefactor). □ 

As a result of ( |2.49| l and Assumption |2.2[ the walk is well-behaved at least along the 
sequence of stopping times Ui, we formalise this in the following result. 

Lemma 2.20. When p is sufficiently close to 1 there exist finite positive constants c and 
C so that for all finite J^-stopping times T with T s {ai : z e N} a.s. and all fc G N 

P(||Xfe-XT|| >w(A:-r)| J-t) a.s. on {T < A:} (2.54) 

and for j < k 

¥{\\Xk-Xj\\> {l + e)s^,^ik-j)\TT) a.s. on {T < j} (2.55) 

with Smax as in Lemma [2.161 


Proof. Note that by Lemma A.l we may assume that T = ai for some I G Z+, the 


general case follows by writing 1 = ^{T=<yi} (a.s.). 

For l|2.54|l, we combine ( |2.49| l and in particular Corollary |2.18| with the proof of 


Lemma 


2.16 For ( |2.55| l, let T' := inf ({ctj : z g M j n [j,oo)) be the time of the next ctj 

shows that P(r' — j > e(fc — j) | Ft) is 


2.19 


after time j. Inequality l |2.52| l from Lemma 
exponentially small in fc — j. On {T' — j < e{k — j)} we use ( |2.54| l starting from time T' 
and simply use the fact that increments are bounded for the initial piece between time 
j and time T'. □ 


For m < rz we say that rz is a (&, s)-cone time point for the decorated path beyond m 


if 


(tube„ U dtube„) n (Z'^ x {—zi, —rz + 1,..., —zzz}) 

C {(x, -j) : zzz < j < rz, ||x - X„\\ < b + s(n- j)}. (2.56) 

In words (see also Figure [^, zz is a cone time point for the decorated path beyond zzz 
if the space-time path together with its i?ioc-tube and decorations by 

determining triangles is contained in cone(6, s, zz — zzz) shifted to the base point (Ar„, —zz); 
recall the definition of cone(b, s, h) in ( |2.25| l. Note that ( |2.56| l in particular implies 

\\Xn —Xj\\<h +s{n —j) for j = zzz,..., zz — 1. (2.57) 


Lemma 2.21. For e > 0, when p is sufficiently close to 1, there exist & > 0 and s > s^ax 
such that for all finite J'-stopping times T with T G {(7^ : z G N} a.s. (i.e., T = aj for a 
suitable random index J) and all A: G N, with T' := inf{cr; : ai> k} 

F(T' is a (6, s)-cone time point for the decorated path beyond T | Ft) > 1 — £ (2.58) 

a.s. on {T < k}. Furthermore 0 < s — Smax 1 can be chosen small. 
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Figure 4: Cone time point n for the decorated path beyond m. Dashed lines indicate the 
i?ioc tube around the path and the triangles are the determining triangles. The cone is 
given by cone(6, s,n — m) shifted to the base point (X„, —n). 


Proof. Denote the event in l |2.57| l (with m = T, n = k and b replaced by 6 — M, where 
M > 0 will be tuned later) by We have 


fc-i 

1 - P(BT.fe I-Ft) < ^P(||Xfe-X,|| >b-M + s{k-j)\j^T) 

j=T 

k—1 k—m—1 

< ^ V{\\Xk-X^\\>b-M\TT)+ V{\\Xk-X^\\>sik-j)\XT). 

j=k—m j—T 

Using ( |2.55| l from Lemma [2 .201 we can make the second sum small by choosing m suf¬ 
ficiently large and s > Smax- Then we can make the first sum small (or even vanish) by 
picking b — M sufficiently large. 

Recall that is the range of the random walk X. Ineguality l |2.52| l from Lemma [2.19| 
implies that P(T' — k > {M — Rioc)/Rk \ Xt) can be made arbitrarily small by choosing 
M sufficiently large. On Bx^k H {T' — k < {M — R\oc)/Rk}, which has high probability 
under P(-1 Xt). we have by construction that 

tuber' n x {-T', -T' + 1,..., -T}) 

C {{x,-j):T<j <T\\\x-XT>\\<b + s{T'-j)}, 


i.e., the path together with its i?ioc-tube is covered by a suitably shifted cone with base 
point {Xt', —T'). 


It remains to verify that under P(-1 Ft) with high probability also the decorations 
(recall ( |2.43| l, ( |2.44| l) are covered by the same cone. T o sho w this we may assume 
T = Ui for notational simplicity; this is justified by Lemma 
as defined in and around l |2.51| l with t = k. 


A.l 


Let Rk 


-(fc) ~(fe) 


be 


Note that dtuber' is contained in a union of space-time rectangles with heights t, 


(fc) 


side lengths 2t^\R\oc V R^) and base points 
A geometric argument shows that on the event 




=(*=)) 


- 




BT,k n Pi > i + 

771=1 


< Ml + eim} U {Rk < i-\- 


m 


})■ 


for d. Ml chosen suitably in relation to b and s, we also have 

dtuber' n x {-T', -T' + 1,..., -T}) 

C {ix,-j):T<j <T\\\x-XT>\\<b + s{T'-j)}. 
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Using inequality l |2.53| l from Lemnia [2.19[ we see that 

oo oo 

P(i?fc >i + m, > Ml + eim | J'r) < XI a.s. on {T < k} 

m—1 m—l 

which can be made arbitrarily small when Mi and ei are suitably tuned. This completes 
the proof of l |2.58| l. □ 

Note that the ai defined in l |2.46| l are themselves not regeneration times since ( |2.49| l 
is in general not an equality of laws. We use another layer in the construction with 
suitably nested cones to forget remaining positive information. 

Recall the definition of cones and cone shells from l |2.25| l, l |2.27| l and FigureThe 
following sets of 'good' w-configurations in conical shells will play a key role in the 
regeneration construction. Let G'(6inn, 6out, Sinn, Sout, h) C {0, ,Si„„, W ,/.) foe the 

set of all w-configurations with the property 

V??o,?7o G {0,lX ^ ^°I-S6out(0) ^ ^ 

UJ G {0,1}^ x{l,...,/t} with a;|cs(&i„„,&out,Sinn,Sout,?i) G G(f'inn, ^out, Sinn, Sout, ^ (2.59) 

r/n(x) = for all (x, n) G cone(6inn, Sinn, h), 

where 77 and rj' are both constructed from l |2.2| l with the same w's. In words, when there 
are I's at the bottom of the outer cone, a configuration from G(6inn, ('out, Smn, Sout, (i) 
guarantees successful coupling inside the inner cone irrespective of what happens out¬ 
side the outer cone. 

Lemma 2.22. For parameters p, 6inn/ ^out, Sinn and Sout as in Lemma [2.13[ 

P(^|cs(feinn,''out,Sinn,Sout,/t) G G((>inn, ('out, Sinn, Sout, ('')) ^ (2.60) 


uniformly in G N. 

Proof. The assertion follows from Lemma [2 .1 3| because if the event GinG 2 defined there 


occurs, then w , 


sih. 


out^Sinn^Sout 


^h) G G(&inn,('out, Sinn, Sout, (i) holds (recall Remark 


2.141. 


□ 


Let us denote the space-time shifts on x Z by i.e., 

0^^,")(t4) = {(x + y, m + n) : (y, m) G A} for A C Z”^ x Z. (2.61) 

An elementary geometric consideration reveals that one can choose a deterministic 
sequence tg 00 with the property that for £ G N and ||x|| < SmaxG+i 

0(°“*^((cone(t^sniax +('out, Sout, f^)) C 0^“ (cOne(5inn, Sinn, G+l)) ■ (2.62) 

Note that this essentially enforces tg « for a suitable p > 1. Indeed, a worst case 
picture (see Figure]^ shows that we need 

f^+is inn “1“ ^inn ~ f^+lS max > tg •^max “I” ^out “I" (tjSout 

which is equivalent to 


ti+i > 


(^(Sout 


-b. 


'^inn '^max 
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Figure 5: Growth condition for the sequence {tg): The small inner cone is cone(t^Si„ax + 
bout, Sout.ti) shifted to the base point (0, —t^). The big outer cone is cone(6inn, Sinn, ft+i) 
shifted to the base point (—f^+iSmax, —The slope of the dashed line is Smax- The 
sequence (f^) must satisfy ux < ci 2 for ux — T ^out T and 0-2 — ^innf^+i T ^inn 

'5 m ax i£+i- 


Thus, we can use (for £ sufficiently large) 

. r r ^ -^out T S max r-n\ 

for any p> -. (2.63) 

^inn ^max 

Furthermore note that using Lemma [2.16| we obtain 

te 

P(3n < te : ||X„|| > Smaxf^) < ^ IP(li-^n|| > Smaxn) < (2.64) 

TL— i Smaxl 

Since te grows exponentially in £, the right hand side is summable in £. Thus, from some 
random £q on, we have sup„<(J|X„|| < s^axte for all £ > £q, and £q has very short tails. 

2.4 Proof of Theorem 12.61 

The proof of Theorem |2.6| relies on a regeneration construction and moment estimates 
for the increments between regeneration times (recall the discussion after Remark [2.7[ l. 
We now have prepared all the ingredients to carry out the argument, which we prefer to 
give in a more verbal, descriptive style. While all the concepts and properties discussed 
below can be easily expressed in mathematical formulas we believe that the resulting 
increase in length and in notational heaviness would burden the text unnecessarily 
without improving neither readability nor understandability. 

That said, the regeneration construction goes as follows (see also Figure [^: 

1. Go to the first ctx after ti, check if 77 in the 6 out-neighbourhood of {X„.,—ai) is = 1, 
the path (together with its tube and decorations) has stayed inside the interior of 
the corresponding conical shell based at the current space-time position and the 
w's in that conical shell are in the good set as defined in l |2.59| l. This has positive 
(in fact, very high) probability (cf. Lemma [2.21[ ) and if it occurs, we have found the 
first regeneration time Tx. 

2. If the event fails, we must try again. We successively check at times £ 2 , t^, etc.: If 
not previously successful, at the f-th step let ae be the first ux after te, check if ae is 
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a cone point for the decorated path beyond with \\X^^ || < Smax^e, the ry’s in the 
6out-neighbourhood of (X^^,—<t£) are = 1, w’s in the corresponding conical shell 
are in the good set as defined in l |2.59| l and the path (with tube and decorations) 
up to time is contained in the box of diameter Soutf^-i + ^out and height 
If this all holds, we have found the first regeneration time Ti. 

(We may assume that is suitably close to this has very high probability 
by Lemma [2 .191 ) 


3. The path containment property holds from some finite £o on. Given the construc¬ 
tion and all the information obtained from it up to the (f—l)-th step, the probability 
that the other requirements occur is uniformly high (for the cone time property 
use Lemma [2.211 with k = ti; use ( |2.49| l to verify that the probability to see rj = 1 
in a box around is high; use Lemma 2.22 to check that conditional on 

the construction so far the probability that the w's in the corresponding conical 
shell are in the good set is high, note that these cu's have not yet been looked at 
by the construction so far). 


4. We will thus at most require a geometric number of tg's to construct the regener¬ 
ation time Ti. Then we shift the space-time origin to {Xti,—Ti) and start afresh, 
noting that by construction, the law of {rj-k-Tiix + 2fTi))xGZ‘^.fcGZ given all the in¬ 
formation obtained in the construction so far equals the law of {Ti-k{x))x&<i,k& 
conditioned on seeing the configuration ryo = 1 in the 6out-box around 0. 

The sequence tg, grows exponentially in I with rate p (see l |2.63| l) and we need to go 
to at most a random t with geometric distribution with a success parameter 1 — <5 
very close to 1. We thus can enforce a finite very high moment of the regeneration 
time; 

P(regeneration after time n) < P(more than log n/log p steps needed) 

1 /I (2.65) 

< jlogn/logp _ ^-a 


where a = log(l/5)/logp can be made large by choosing 6 small and p close to 1. 
Both is achieved by choosing p close to 1. 


We obtain a sequence of random times Ti < T 2 < ■ ■ ■ such that (Xri — — 

7i-i)i=2.3,... are i.i.d. and E[T^^],E[(T 2 — Ti)*'] < 00 and hence also E[||2fTi||^],E[|| 2 fT 2 — 
< 00 for some b > 2. The existence of such regeneration times implies The¬ 
orem [2]^ by standard arguments, see e.g. the proof of Corollary 1 in Il20l (it is easy to 
see from the construction that Xt^ — is not a deterministic multiple of Ti — Ti_i) 
and the proof of Theorem 4.1. in II28II for the functional CLT. Note that the speed must 
be 0 by the assumed symmetry; see Assumption |2.3[ 


Remark 2.23. In the general case without the Assumption 2.3 the above argument 
yields that there must be a limit speed, its value is given only implicitly as E[ArT 2 — 
a:tJ/e[T 2 - Ti]. 

If in Assumption |2.3| we would additionally require symmetries with respect to co¬ 
ordinate permutations and with respect to reflections along coordinate hyperplanes 
then the limiting law $ would be a (non-trivial) centred isotropic d-dimensional normal 
law, cf. the proof of Theorem 1.1 in |[3. 


3 A more abstract set-up 

The goal of this section is to present an abstract set-up where a renewal construction 
similar to the one of the previous section can be implemented. Our main motivation of 
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Figure 6 ; A schematic example; The walk passing through a sequence of cones in an 
attempt to regenerate. Here, n = tq + fa- 


this set-up is to study the dynamics of ancestral lineages in spatial populations, but it 
can be likely applied for other types of directed random walks in random environment. 

In Sections |3.1| and |3.2[ we present certain abstract assumptions on the random 
environment and the associated random walk. These assumptions allow to control the 
behaviour of the random walk using a regeneration construction that is very similar 
to the one from Section In particular, they allow to link the model with oriented 
percolation, using a coarse-graining technique. 

We would like to stress that coarse-graining does not convert the presented model 
to the one of the previous section. In particular, the nature of regenerations is some¬ 
what different. We will see that the sequence of regeneration times and associated 
displacements, (T^+i — — XTi)i >2 is not i.i.d. but can be generated as a certain 

function of an irreducible, finite-state Markov chain and additional randomness. By 
ergodic properties of such chains, this will lead to the same results as previously. 

Theorem 3.1. Let the random environment 77 and the random walk X satisfy the as¬ 
sumptions of Sections |3 .1 1 and |3 .2 [ below with sufficiently small parameter eu- Then the 
random walk X satisfies the strong law of large numbers with speed 0 and the annealed 
central limit theorem with non-trivial covariance matrix. A corresponding functional 
central limit theorem holds as well. 

A concrete example satisfying the abstract assumptions of Sections |3.1| and |3.2| will 
be given in Section They can also be verified for the oriented random walk on the 
backbone of the oriented percolation cluster which was treated in O using simpler, but 
related, methods. 

3.1 Assumptions for the environment 

We now formulate two assumptions on the random environment. The first assump¬ 
tion requires that the environment is Markovian (in the positive time direction), and 
that there is a 'flow construction' for this Markov process, coupling the processes with 
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different starting conditions. The second assumption then allows to use the coarse- 
graining techniques and the links with oriented percolation. 

Formally, let 

U '.= {U{x,n) : X S n G Z} 

be an i.i.d. random field, C/(0,0) taking values in some Polish space Eu {Eu could be 
{—1, +1}, [0,1], a path space, etc.). Furthermore for G N let (0) c Z'^ be 

the ball of radius i?,, around 0 with respect to sup-norm. Let 

ip : Z^”” X Z+ 

be a measurable function. 

Assumption 3.2 (Markovian, local dynamics, flow construction). We assume that 77 := 
(r 7 n)nGZ is a Markov chain with values in Z^ whose evolution is local in the sense that 
? 7 „+i(x) depends only on r]n{y) for 1 / in a finite ball around x. In particular we assume 
that 77 can be realised using the 'driving noise’ U as 

77 „+i(x) = (/?( 6 »'" 77 „|^^^ , 6I^(7( • ,n + 1)1^^^ ), x G Z"*, G Z. (3.1) 

Here 9^ denotes the spatial shift by x, i.e., 9^r]n{ ■) = 77 „( • + x) and 6 ^U{ • , 7 r + 1) = 
U{- +x,ri + l). Furthermore 9^r]n\Bii + ^)\br are the corresponding 

restrictions to the ball Br^ ■ 

Note that ( |3.1| l defines a flow, in the sense that given a realisation of U we can 
construct 77 simultaneously for all starting configurations. In most situations we have in 
mind the constant zero configuration 0 G Z^ is an equilibrium for 77 , that is, 

’ ■) = 

and there is another non-trivial equilibrium. It will be a consequence of our assumptions 
that the latter is in fact the unique non-trivial ergodic equilibrium. 

The second assumption, inspired by ®, allows for comparison of 77 with a supercrit¬ 
ical oriented percolation on a suitable space-time grid. Loosely speaking, this assump¬ 
tion states that if we have a good configuration on the bottom of a (suitably big) block 
and the driving noise inside the blocks is good, too, then the configuration on the top of 
the block is also good and the good region grows with high probability. Furthermore if 
we input two good configurations at the bottom of the block then good noise inside the 
block produces a coupled region at the top of the block. 

Formally, let Lt,Ls G N. We use space-time boxes whose 'bottom parts' are centred 
at points in the coarse-grained grid LsZ‘^ x LtZ. They will be partly overlapping in the 
spatial direction but not in the temporal direction, and we typically think of Lt > Lg ^ 
R,,. 

For (x, n) G Z'^ X Z we set 

blockm(x, n) := {(j/, fc) G Z'^ X Z : \\y — Lsx\\ < mLs,nLt < k < {n + (3.2) 

and block(x, n) := blocki(x, ri); see Figure For a set zl c Z'^ x Z, slightly abusing 
the notation, we denote by U\^ the restriction of the random field U to A. In particu¬ 
lar, ff|biock 4 ( 2 n) the restriction of U to block 4 (x, n) and can be viewed as element of 

Eu 

Assumption 3.3 ('Good' noise configurations and propagation of coupling). There ex¬ 
ist a finite set of 'good' local configurations c and a set of 'good' local 

realisations of the driving noise Gu C ^-j^g following properties: 
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(n + l)Lt 


nLt 


(x Kr{)Ls 


xLs 


{x + Krf)La 


Figure 7: Locality of the construction of (? 7 „) on the block level for d = 1. If 17 is known 
in the grey region and rjnL^ is known on the bottom of the dashed trapezium then the 
configurations 77 ^ are completely determined inside block(x, n) drawn in solid lines. 


• For a suitably small eu. 


^Iblock4(0,0) S ^u) >!-£[/ 


(3.3) 


• For any {x, n) x Z and any configurations rinLt,v'nLt ^ 

’ ^nLt\ B2La(Lax) ^ ^lblock4(£,n) ^ 

=> V(n+i)L,{y) = v[n+i)LSy) for all 7/with ||?/- Ls^|| < 3Ls (3-4) 

and V(n+i)L,\B,^^^Laix+^) ^ for all e with ||e|| < 1 , 


and 

?7KLtlB24,^(L,5r) = ynulB^Lai^sx) ^ = ^kiv) ^r all {y,k) e block(J,n), 

(3.5) 

where 77 = ( 77 ^) and 77 ' = ( 77 ^) are given by ( |3.1| l with the same U but possibly 
different initial conditions. 

• There is a fixed (e.g., Lg-periodic or even constant in space) reference configura¬ 
tion 77 ™^ G Z^ such that x) ^ x &Z'^. 

Note that if the event in l |3.4| l-l [331 l occurs then a coupling of 77 and 77 ' on i? 2 L^ {Lsx) x 
{nLt} has propagated to L2Ls(Ls(2:+e))x {(n+l)Lt} for l|e|| < 1 and also the fact that the 
local configuration is 'good' has propagated. The event in ( |3.4| l enforces propagation 
of goodness and can also be viewed as a contractivity property of the local dynamics. 
In other words the flow tends to merge local configurations once they are in the 'good 
set'. 

Remark 3.4 (Locality on the block level). Put 

:=L^(r^l+1). (3-6) 

From the local construction of 77 given in l |3.1| l it follows easily (see Figure [7]l that for 
fixed (x, n) G x Z the values 77„(x) for (x,n) G block(x, n) are completely determined 
by 77 KLt restricted to Bk^Ls{xLs) and U restricted to U||j 7 ||<if^block(x + y,n). 
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Using the above assumptions, it is fairly standard to couple 77 to an oriented percol¬ 
ation cluster. Recall the notation in Section [2.1 [ and in particular the definition of the 
stationary discrete time contact process in l |2.3| l. 

Lemma 3.5 (Coupling with oriented percolation). Put 

= 0 - 7 ) 

If eu is sufficiently small, we can couple U(x, n) to an i.i.d. Bernoulli random field u(x, n) 
with ¥{Lj(x,n) = 1 ) > 1 — cs such that (7 > w, and ej can be chosen small (how small 
depends on eu, of course). 

Moreover, the process 77 then has a unigue non-trivial ergodic equilibrium and one 
can couple a stationary process 77 = {r]n)nez with 770 distributed according to that equi¬ 
librium with w so that 


G{x, n) := U{x^ 7 z)l 






eG,,} 


>^{x,n), (a;, 77 .) C X Z 


(3.8) 


where ^ := {^(a;, n) : x G Z*^, n c Z} is the discrete time contact process defined by 

^{x,n) ■■= l{zdx{-oo}-)-"(a;,n)}' 

Proof. The first part is standard; Note that the U(x,n)'s are i.i.d. in the n-coordinate, 
with finite range dependence in the ^-coordinate. Using l |3.3| l, l |3.4| l and ( |3.5| l, we can 
employ e.g. the Liggett-Schonman-Stacey device f ll2^ or ll^ Thm. B26]). 

For the second part consider for each k G N the process = {vin'^)n>-kLt which 
starts from and evolves according to for n > —kLt, using given U's 

which are coupled to uj's as above s o tha t U > w ho lds. We see from the coupling 
properties guaranteed by Assumption 3.3 and Lemma 3.12 below that the law of 
restricted to any finite space-time window converges. By a diagonal argument we can 
take a subsequence km Z' 00 such that 77 „(a:) := lim^-^oo exists a.s. for all 

(a;, n) s Z'^ X Z, then l |3.8| l and ( |3.9| l hold by construction. 

The fact that the law of limit is the unique non-trivial ergodic equilibrium can be 
proved analogously to JH Cor. 4]. □ 


Remark 3.6 (Clarification about the relation between ^ and 77 ). The contact process ^ is 
defined here with respect to w analogously to the definition of the discrete time contact 
process 77 with respect to w in ( |2.3| l. The rationale behind this change of notation is that 
throughout the paper 77 is a stationary population process (contact process in Section 
and logistic BRW in Section and the random walk X is interpreted as an ancestral 
lineage of an individual from that population. The coarse-grained contact process ^ 
plays a different role. In particular, the knowledge of ^ alone does not determine the 
dynamics of X; cf. definition of X in ( |3.11D . 

Finally, we need the following technical assumption which is sufficiently strong for 
our purposes but can be relaxed presumably. 

Assumption 3.7 (Irreducibillty on G^). On G{x,n), conditioned on seeing a particular 
local configuration x G Z^^^= n Gjj at the bottom of the space-time box [time co¬ 
ordinate nLt], every configuration x' G Gn has a uniformly positive chance of appearing 
at the top of the space-time box [time coordinate (n -I- l)Lt]- 

Remark 3.8. For the discrete time contact process the above assumptions can be 
checked easily in the case d = 1 when p is sufficiently close to 1. For G^ we could 
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for instance take configurations rj G { 0 , 1 }^”^ with 


#{||a:|| < Ls/2,r]{x) = 1} > -L^. 


For Gu we could take configurations of w’s for which this property propagates to the 
top of the block and its neighbours irrespective of the positions of the I's in the initial 


configuration (cf. construction in the proof of Lemma 2.131. For d > 2 one can reduce 
the argument to the one-dimensional case. 


3.2 Assumptions for random walk 

We now state the assumptions for the random walk X = {Xk)k=o,i,... in the random 
environment generated by 77 . To this end let U := {U{x,k) : x G IJ^,k G Z_|_) be an 
independent space-time i.i.d. field of random variables uniformly distributed on ( 0 , 1 ). 
Furthermore let 

X X [0,1] ^ Br^ (3.10) 

a measurable function, where Rx G N is an upper bound on the jump size as well as on 
the dependence range. Given 77, let Xq = 0 and put 

Xk+i.= Xk + ^x{e^^r^-k\j,^^,e^>‘r^_k-i\g^^.U{Xk,k)), fc = 0,l,.... (3.11) 

Note that, as usual, forwards time direction for X is backwards time direction for 77. 
Assumption 3.9 (Closeness to SRW while on G = 1). A walker with dynamics ( |3.11| l 
starting from the middle half of the top of a box with G{x, n) = 1 stays inside the box 
with high probability: 

min f( max W^k — 2 :|| < —p 
z •. \\z—x\\<Ls/2 V (n—l)Lt<fc<nLt 4 

(3.12) 


X, 


(n-l)Lt 


= z,G(a;,n) = l, 77 j 


> 1 - e. 


Remark 3.10. (a) Note that ( |3.12| l translates intojhe upper bound eRx + Ls/(4Lt) for 
the speed of the walk X on a block satisfying G{x, n) = 1. The factor j in ^ inside 
( |3.12| l is somewhat arbitrary. Depending on eu and the ratio of Lg to Lt one could 
use a different factor. 


(b) The simple Assumption |3 . 9| allows to obtain a rough a priori bound on the speed of 
the walk and suffices for our purposes here, a more elaborate version would reguire 
successful couplings of the coordinates of X with true random walks with a small 
drift while on the box, similar to the proof of Lemma [2.16[ 

Assumption 3.11 (Symmetry of ipx w-f-t. point reflection). Let g be the (spatial) point 
reflection operator acting on 77, i.e., Qrjk{x) — r]k{—x) for any fc G Z and x G We 
assume 

^x{grio\Bj^^ > ,^( 0 , 0 )) = ( 3 . 13 ) 

Note that l |3.13| l guarantees that the averaged speed of X will be 0. 


3.3 The determining cluster of a block 

We now explain how Theorem |3.1| can be proved using similar ideas as in Section 
In order to avoid repetitions and to keep the length of the paper acceptable, we only 
explain the major differences to the proof of Theorem|2.6[ 
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The main change that should be dealt with is the fact that the construction of the 
random walk X requires not only the knowledge of the coarse-grained oriented percol¬ 
ation but also of the underlying random environment 77 . This additional dependence 
on rj should be controlled at regeneration times. To tackle this problem, Assumption |3.3| 
and Lemma [3^ play the key role. By this lemma, the value of r]{x, n) can be reconstruc¬ 
ted by looking only at the driving noise U in certain finite set 'below' (a;, n). 

Formally, for (a?, n) s x Z we define its determining cluster DC (a?, n) by the follow¬ 
ing recursive algorithm: 

1. Initially, put k := n, DC(a:,n) := {(a:,n)}. 

2. If ^(y, k) = \ for all (y, k) € DC(a;, n) : Stop. 

3. Otherwise, for all blocks where this condition fails, add every block one time layer 
below that could have influenced it (cf Remark [3.4[ l, that is replace DC(a?, n) by 

DC(a;, n) U {(z, fc — 1) : \\z — y|j < for some y with ^(y, k) = O}, (3.14) 

put k ■— k — I and go back to Step 2. 

Lemma 3.12. For eu small enough, the height (and the diameter) of DC(a:, n), defined 
as 

height(DC(a;, n)) := max {n — k : (y, k) S DC(a:, n)}, (3.15) 

is finite a.s. with exponential tail bounds. 

Proof. This can be shown as in the proof of Lemma |2.11| see alternatively Lemma 7 
in llllll . or proof of Lemma 14 in |[ 6 l. □ 

Remark 3.13. On{^(ir,n) = 1}, r]\biock(xn) is a function of local randomness. In fact it is 
then determined by (^lbiock 5 (s,K)ubiock 5 ( 5 :,f;-i)- Thus, 77 on block(x,n) is determined by the 
'wet boundary' plus local randomness in a slightly 'thickened' version of DC{x,n) which 
is the analogue of the 'determining triangle' D{x,n) from ( |2.43| l in this coarse-grained 
context. 

To see this, consider the system 77 ' := ( 77 ^ : (n—l)Lt <n< {n+l)Lt) which starts from 
v[n-i)Lt — tbe fixed boundary condition 77 ^(y) = 77 ''®^(y) for ||y — Lsx\\ > 5Ls 

and (n — l)Lt < n < (n + l)Lt- For (y, n) S block 5 (a;, n) U block 5 (a:, n — 1) the values 77 ^(y) 
are computed using ( |3.1| l with the same realisations of U as the true system 77 . 

Note that ^{x,n) = 1 implies that S Gu and 

U\b\ocki{x+e,n-l))^Gr,{ | (L,, (2+5)) ) ~ ^ 

Now use l |3.4| l to see that and -q^Lt agree on B 2 lX^sx), then use that and ( |3.4| )-( [T5l l 
to verify that 77 ' and 77 agree on block(x, n). 

3.4 A regeneration structure 

In this section we construct regeneration times similar to those constructed in Sec¬ 
tion [2]^ First we need to introduce the analogue of the 'tube around the path' and its 
'decoration with determining triangles'; cf. equations l |2.42| l, ( |2.43| l and ( |2.44| l. We set 



:= {x : 3k, {m - l)Lt <k < mLt, \\Xk - xLt\\ < Ls + Rx}, 

(3.16) 

Tube^ 

■= IJ T* X {to}, 

(3.17) 


rnKn 


DTube^ 

:= U DC(2?,j). 

(3.18) 


(5,J)GTube~ 
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We define the coarse-graining function tt : by 


'k{x) = 7r(xi, ■■.,Xd) 




(3.19) 


and denote by p{x) the relative position of x inside the block centred at xL^, i.e. we set 


p{x) := X — xLs- 


(3.20) 


We define the coarse-grained random walk X = {Xn)n=o,i,... and the relative positions 

y = (yn)n=o.,i.... by 


X^:=li(X7iLj and := p(^nLj. (3.21) 

We need to keep track of the relative positions to preserve the Markovian structure. 
Note that between the original random walk and the coarse-grained components just 
defined we have the following relation: 

^nLt ~ Xn' 


We define the filtration T := (Jk)K=o,i,... by 
Jk := (7(^{Xj,Yj) : 0 < J< n) V cr(w(y, j),j), C^lblock 4 (y, j) ' ^ DTubejiy (3.22) 

To mimic the proofs of Section for the model considered here we need the following 
ingredients: 


1. As in Lemma 


2.16 


there exist Smax (that is close to | under our assumptions) and 


positive constants C, c such that 

P(l!^n|| > wn) < Ce-™. 

2. For stopping times (analogous to a's in ( |2.46| l) we set 

Dn :=n-\- max |height(DC(2:, n)) : x G Vk| 

and define 


lin I 


(To ■= 0; cTi := min \ m > di-i : _ max _ < m [•, i > 1. 

cri_i<n<m 




(3.23) 


(3.24) 


(3.25) 


Lemma 3.14. When 1 — is sufficiently close to 1 there exist finite positive constants 
c and C so that 


P((Ti+i — ai> n \ Xsi) < Ce for all n = 1,2,..., i = 0,1,... a.s., (3.26) 

in particular, all Oi are a.s. finite. Furthermore, 

^((w(-,-j-Ui)j=o,i,... ^ for every z = 0,1,... a.s., (3.27) 

where V' denotes stochastic domination. 

Proof. Analogous to the proof of Lemma [2.171 (see also Lemma[37T^. □ 
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Similarly to the definition in ( |2.56| l we say that n is a {b, s)-cone time point for the 
decorated path beyond m (with m < n) if 

DTubes n x {—n, —n + 1,..., —m}) 

r ~ ^ ^ ^ J ~ ^ , (3.28) 

c {(x,-j) : TO < j < n, ||a; - XkII < 6 + s(n - j)}. 

In words, as in Section[2]^n is a cone time point for the decorated path beyond to if the 
space-time path {Xj, together with its 'tilde'-decorations is contained in the 

cone with base radius b, slope s and base point —n). 

Lemma 3.15. There exist suitable b and s > S^max such that for all finite J^-stopping 
times T with T G {ai : f G N} a.s. (i.e., T = aj for a suitable random index J) and all 
fc G N, with T' := inf{(Ti : Ui > fc} 

P(r' is a (5, s)-cone time point for the decorated path beyond T | Xf) > 1 — c (3.29) 

a.s. on {T < k}. Furthermore 0 < s — s^^ax ‘C 1 can be chosen small. 

Proof. Analogous to the proof of Lemma |2.21[ Intermediate results, that is Lemma [2.19| 
and Lemma [2.201 can be adapted to the present situation. □ 

We now define 'good configurations' of w's (analogous to l |2.59| l). Recall the defini¬ 
tion of a cone shell in ( |2.27| l. Let G'(&inn, bout, Sinn, Souuh) C {Q, be 

the set of possible w-configurations in cs(&inn, bout, Sinn, Sout, h) with the property 

€ {0,1}^" = 1 and 

^ ^ {^1 1} ^ ^ ‘^|c5(6inii,^out,Sinn,Sout,^) ^ G^(^inn; ^OUt 5 5111115 ^OUt 5 ■ (3.30) 

^(x, n) = ^'{x, n) for all {x, n) G cone(5in„, Sinn, h) 

where ^ and are both constructed from time 0 using the same w's, i.e. when A and A' 
are subsets of with = ^(-,0) and 1^' = C^(-,0) then (cf. ( |2.1| l) 

Note that if f,{x, 0) = 1 in the ball (0) and 

^|cs(binn,&out,Sinn,Sout A) ^ i ^^OUt ; ^inil; ^OUt j 

then 


{r]n{x) : {x,n) G block(a;,n), (x,n) G cone(6inn, Sinn, ^)} 


is a function of 770(2/), ||j/|| < ^out^s and ( 7 |biock 4 (s,n)' ^ cone( 6 inn,Sinn, h). In particu¬ 

lar, if we start with different t/q and U' with Voiv) = Vo{y), \\y\\ < (^out^s and C/'lbiock4(5.K) = 

^\h\oc'ki{x,n)' (^i ^ COne(5inni Sinni then 

ynix) = r]'„{x) for all {x,n) G block(x, n), (a?, n) G cone(&inn, Sinn, ^)- 


Proof sketch for Theorem 13. II We now have all the ingredients for the regeneration 
construction, to imitate the proof of Theorem |2.6[ We again choose to keep the argu¬ 
ments more verbal and descriptive, hoping to strike a sensible balance between nota- 
tional precision and readability. 


First we choose a seguence to,ti,... with f 00 such that ( |2.62| l is satisfied with 
ax replacing s 


and p arameters 6out. Sout. (>inn and 


Recall from Remark 

(^1 blocks (5 ;,t/)U blocks (x,n—l) ' 


3.13 


that on the event {C(a;,n) = 1}, ?7lbiock(: 


adapted from Lemma 3.15 
^ is determined by 
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1. Go to the first ai after ti, check if in the 6out-neighbourhood of we 

have ^ = 1, the path (together with its tube and decorations) has stayed inside 
the interior of the corresponding conical shell based at the current space-time 
position and the w's in that conical shell are in the good set as defined in ( |3.30| l. 
This has positive (in fact, very high) probability (cf. Lemma [3.151 1 and if it occurs, 
we have found the 'regeneration time'. 


2. If the event fails, we must try again. We successively check at times t 2 , etc.: If 
not previously successful, at the £-th step let aj(e) be the first after ti, check if 
is a cone point for the decorated path beyond ti-i with |j < SmaxCTj(£), 

the ? 7 ’s in the 6out-neighbourhood of are = 1, w's in the corresponding 

conical shell are in the good set as defined in ( |3.30| l and the path (with tube and 
decorations) up to time is contained in the box of diameter Soutf^-i + &out and 
height f^_i. If this all holds, we have found the regeneration time. 

(We may assume that is suitably close to ti-i, this has very high probability 

by an adaptation of Lemma [2TT^) 


3. The path containment property holds from some finite on. Given the construc¬ 
tion and all the information obtained from it up to the l)-th step, the probability 
that the ot her re quirements occur is uniformly high: For the cone time property 
with k = tr, use ( |3.27| l to verify that the probability to see ^ = 1 in 


use Lemma 


3.15 


2.22 


a box around is high; use (a notational adaptation of) Lemma 

to check that conditional on the construction so far the probability that the w’s in 
the corresponding conical shell are in the good set G(6inn, ^out, Smn, Sout, te) is high. 
Note that these w's have not yet been looked at. 


4. We thus construct a random time i?i with the following properties: 


(ii) the decorated path up to time Ri is in cone(6inn, Smn, Ri) centred at , Ri); 

(hi) after centring the cone at base point b t s s t Ri) 

the good set ^out? .^l)- 

We will thus at most require a geometric number of t^'s to construct the i?i. As in 
step 1^ in the proof of Theorem |2.6| we obtain 

P(i?i > n) < P(more than log n/log c steps needed) < ^'°s"/iogc _ ~-a^ 


where again a can be chosen large when p is close to 1. 
5. Set 


Vl ■— ■ ll^ll — ^outl^s)) 

Yi := Yg^,the displacement of X^^^^ relative to the centre 
of the Lg-box in which it is contained. 

Now we shift the space-time origin to {Xj^^,Ri) (on coarse-grained level). Then 
we start afresh conditioned on seeing 

(i) configuration ^ = 1 in the &out-box around 0 (on the coarse-grained level); 

(ii) 771 on the bontLs box (on the 'fine' level); 
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(iii) Displacement of the walker on the fine level relative to the centre of the 
corresponding coarse-graining box given by Yi. 

6. We iterate the above construction to obtain a sequence of random times Ri, posi¬ 
tions , relative displacements Y and local configurations By construction 

is a Markov chain. Furthermore, {Yi,rji)^^^ is itself a finite state space Markov 
chain and the increments — Xj^ ,Ri+i — Ri) depend only on {Yi,%). 

Along the random times L^Rn, 

n 

^LtR„ = ~ ^Ri-^ 

i=l 

is an additive functional of a well-behaved Markov chain (with exponential mixing prop¬ 
erties) and 


E[{R,+,-R,r\%rji\ <oo, < cx) 

for some a > 2 uniformly in Yi,% (cf. Step 4). From this representation the (functional) 
central limit theorem can be deduced; see e.g. Chapter 1 in IfT^ or Theorem 2 in 026II . 

Note that the speed of the random walk must be 0 by the symmetry assumption; see 
( [3T3l l. □ 

4 Example: an ancestral lineage of logistic branching 
random walks 

In this section we consider a concrete stochastic model for a locally regulated, spatially 
distributed population that was introduced and studied in lO and we refer the reader 
to that paper for a more detailed description, interpretation, context and properties. 
We call this logistic branching random walk because the function / in l |4.1| l, which 
describes the dynamics of the local mean offspring numbers, is a 'spatial relative' of the 
classical logistic function x x{\ — x) which appears in many (deterministic) models for 
population growth under limited resources. See also Remark below for a discussion 
of related models and possible extensions. 

After defining the model we recall and slightly improve some relevant results from 
Il6ll . Then in Proposition |4.7| we show that in a high-density regime (see Assumption |4.2| l 
assumptions from Section are fulfilled by the logistic branching random walk and the 
corresponding ancestral random walk. 

4.1 Ancestral lineages in a locally regulated model 

Letp = iPxy)x,yGZ'^ = iPy-x)x,yGZ‘^ ^0 a Symmetric aperiodic stochastic kernel with finite 
range Rp > 1. Furthermore let A = {Xxy)x,yeZ‘^ be a non-negative symmetric kernel 
satisfying 0 < Xxy = Ao,y-x and having finite range Rx- We set Aq := Aqo and for a 
configuration ( e and x we define 

fix; C) := Cix){m - AoC(a:) - XxzCiz))^ ■ (4-1) 

Z^X 
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We consider a population process 77 := {r]n)„ez with values in , where as in the 
previous sections rin{x) is the number of individuals at time n G Z at site x G Z^. 
Before giving a formal definition of ry let us describe the dynamics informally: Given the 
configuration in generation n, each individual at x (if any at all present) has a Poisson 
distributed number of offspring with mean f{x;r]n)/r]n{x), independent of everything 
else. Offspring then take an independent random walk step according to the kernel p 
from the location of their mother. Then the offspring of all individuals together form 
the next generations configuration 7 y„+i. For obvious reasons p and A are referred to as 
migration and competition kernels respectively. Note that in the case A = 0 the process 
r] is literally a branching random walk. 

We now give a formal construction of ry. Let 

U := :nGZ, x,yGZ‘^, ||a: - y|| < Rp} (4.2) 

be a collection of independent Poisson processes on [0, 00 ) with intensity measures of 
given by py^ dt. The natural state space for each Un'^'’ is 

I? := {'0 : [0, 00 ) —>■ Z+: y/i cadlag, piece-wise constant only jumps of size l}, (4.3) 

which is a Polish space as a closed subset of the (usual) Skorokhod space V. For given 
Vn S Z^, define ry„ +1 G via 

Pn+i{x)-.= Y. uiy’^\f{y.Vn)), XGZ'^. (4.4) 

y '■ \\x-v\\<Rp 

Note that for each x, the right-hand side of l |4-4| l is a finite sum of (conditionally) Poisson 
random variables with finite means bounded by ||/||oo- Thus, ( |4.4| l is well defined for 
any initial condition - in this discrete time scenario, no growth condition at infinity, etc. 
is necessary. Furthermore we note that by well known properties of Poisson processes 
7y„+i, given 7y„, is a family of conditionally independent random variables with 

ry„+i(a;) ~ Pois(^^ Py 2 ;/(y;ry„)^, xGZ‘^. (4.5) 

yCJA 


For —00 < m < n set 

Gm,n '■= : m < k < n, x,y G Z‘^). (4.6) 

By iterating l |4.4| l, we can define a random family of ^^.n-measurable mappings 

-cx) <m <n such that 7y„ = (4.7) 

To this end define ‘f’m.m+i as in l |4.4| l via 

($^,™+i(C))(x) := ^ Ut-\fiy;0) for y e and C G (4.8) 

V ■ \\x-y\\<Rp 

and then put 


4^m,n — l,n O • • • O 4*^ (4.9) 

Using these mappings we can define the dynamics of {r]n)n=m,m+i,... simultaneously for 
all initial conditions rym G Z^ for any to G Z. 

Let us for a moment consider the process ry = (ryn)„=o,i,...- Obviously, the configura¬ 
tion 0 G Z^ is an absorbing state for ?y. Thus, the Dirac measure in this configuration is 
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a trivial invariant distribution of rj. In (61 it is shown that for certain parameter regions, 
in particular m G (1,4) and suitable A, the population survives with positive probabil¬ 
ity. For m G (1,3) (and again suitable A) the existence and unigueness of non-trivial 
invariant distribution is proven. We recall the relevant results for m G (1,3). 

Proposition 4.1 (Survival and complete convergence, O). Assume m G (1,3) and letp 
and A be as above. 

(i) There are Ag = Ao(m,p) > 0 and a* = a*{m,p) > 0 such that if Ag < Ag and 

< a*Xo then the process {r]n)n=o,i,... survives with positive probability (if 
survival for one step has positive probability) and has a unigue non-trivial invariant 
extremal distribution P. 

(ii) Conditioned on non-extinction, rjn converges in distribution in the vague topology 
to P. 


Since we are only interested in the regime when the corresponding deterministic 
system, cf. I |4.14| l below, is well controlled and in particular. Proposition |4. 1 1 guarantees 
that a non-trivial invariant extremal distribution P exists, we make the following general 
assumption. 


Assumption 4.2. 1. With the notation from Proposition 4.1 

and X; x^O ^ U Ag. 


we assume m 


G (1,3) 


2 . 7 := J2x ^ox is sufficiently small. 

Note that a* is determined by the dimension d, the parameters to, p and a renorm¬ 
alised A by the reguirement that the left-hand side of ( |4.19| l at (^ = m* must be strictly 
smaller than 1, see Section [4?^ below. 

Under this assumption we can (and do so from now on) consider the stationary 
process 77 = {r]n)n& with ? 7 „ distributed according to P. From the informal description 
(after ( |4.1| l) above and the formal definition ( |4.4| l it is clear that the model can be easily 
enriched with genealogical information; see e.g. Chapter 4 in O. Put 


Pj^{k;x,y) 


Pyxfiy;V-k-i) 

Y.zPzxf{z-,V-k-i)’ 


x,y k € Z+ 


(4.10) 


with some arbitrary convention if the denominator is 0. For a given 77 , conditioned on 
77 g( 0 ) > 0, let X := {Xk)k=o,i, 2 ,... be a time-inhomogeneous Markov chain with 


Xg = 0, and F{Xk+i = y\Xk = x,r]) = Pr,{k-,x,y). (4.11) 


This is the dynamics of the space-time embedding of the ancestral lineage of an in¬ 
dividual sampled at random from the (space-time) origin at stationarity, conditioned 
on the (full) space-time configuration 77 . Note that given 77 , we see from l |4.4| l that 
the number of offspring coming from y in generation —k — 1 that moved to x is given 
by {f{y;ri-k-i)) which is Pois(pya;/(7/; 77_fe_i))-distributed conditional on the sum 
over all y in the neighbourhood of x being equal to ri-k{x). Since a vector of independent 
Poisson random variables conditioned on its total sum has a multinomial distribution we 
see that the dynamics of the ancestral lineage are indeed given by ( |4.11| l. 

Our main result in this section is the following theorem. 

Theorem 4.3 (LLN and averaged CLT). Assume d > 1, let the Assumption |4.2 | be satisfied 
and let 77 = {ri„)nez be the stationary process conditioned on 77 g( 0 ) > 0. For the random 
walk (Alfc)fe=g.i... defined in ( |4.11| l we have 

Py(^^Xk ^ 0 '^ = 1 for P(-|77 o(0) >0)-a.a. 77 , (4.12) 
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and for any g e C';,(K‘^) 


E 


'g{Xk/Vk) |r?o(0) >0] 


(4.13) 


where $ is a non-trivial d-dimensional normal law and <i)(g) := J g{x) 4>(dx). 


Proof. The assertions of the theorem follow from a combination of Proposition |4.7| and 
Theorem 13. II □ 


4.2 Deterministic dynamics 

For comparison, we consider the dynamical system (such systems are also called a 
coupled map lattices) C := (Cn)n=o.i,... on [0,oo)^ defined by 

Qn{x) ^ Pyxf{yXu-i), xSZ"', neN (4.14) 

with / from and arbitrary initial condition Co € [0, (cf. 0 Eq. (5)]). It is easily 
seen from ( |4.1| l that with 

777 — 1 

= , (4.15) 

C*(-) = m* is an equilibrium of the dynamical system C- Furthermore, setting 

T •” ^ ^ ^xy ■— ^xy!^ and Cn(^) ■— 'yCn{x') 

z 

we see from (|43 that (Cn)n=o,i.... solves 

Cn{x) = ^PyxCn-liy){rn-XoCn-liy)-^\zCn-liz))~'', (4.16) 

j/GZ'i z/a: 

i.e., ( I4.14D with A in the function / replaced by A. Thus, we can and shall assume 7 = 1 
for the rest of this subsection. 

Lemma 4.4. There exist oq < a < m* < fi, e = e{m, A) > 0, Rq, ko, Nq and sq such that 
for all R> Ro the following assertions hold: 

(i) If Co(y) G [a,/3] for all y e Br{x) then 

Cn(2/) e [(1 +e)a,/3/(l +e)] for all n > Nq, \\y - x\\ < R + so{n- Nq), (4.17) 

and 


Cn{y) > Oo for all u > 1 , \\y - x\\ < R - ko + squ. 


(ii) For iCiy))yeBu^{x) G [a ,we have 


E 

{x) 


d 


f{x;0 



(4.18) 


(4.19) 


(hi) One can choose P — a > 0 arbitrarily close to 0. 

Proof. Assertions (i) and (hi) follow from Lemma 11 and 12 (and arguments in their 
proofs) in O. For (ii) see the proof of Lemma 13 and in particular Eq. (40) in O. □ 

Remark 4.5 (Interpretation of LemmajA^. Assertion (i) in the above lemma means that 
if Co in the neighbourhood of x is in the interval [a, P] around m* then the regions around 
X where Cn is bounded away from 0 and where it is close to m* grow at positive speed 
(after a finite number of steps). Assertion (ii) means that the equilibrium C*(-) = ra* is 
attracting. 
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4.3 Coupling reloaded 


Remark 4.6 (Initial/boundary cond itions on certain space-time regions). Note that for 


any n G N, $o,n as defined in (4.8 1 can be viewed as a function of {U, 

n, x,y G 

Let L g'H, Rp the range of p, put 


{x,v) 


: 0 < m < 


cone 


(L, Rp) := {(a;, n) G x : ||a;|| < L -I- i?pn} 


(4.20) 


(recalling ( |2.25| l, we have cone(L, i?p) = U/i>ocone(L, i?p, h)). For given values of r]kix), 
{x,k) G X Z+) \ cone(L,iip)) U ([—L, x {0}) (we can view the latter set as a 
'space-time boundary' of cone(L, Rp)), we can define consistently inside cone(L, Rp) 
through ( |4.4| l. 

In fact, we can think of constructing the space-time field ry in a two-step procedure: 
First, generate the values outside cone(L, Rp) (in any way consistent with the model), 
then, conditionally on their outcome, use ( |4.4| l inside. 


Proposition 4.7. Let Assumption 4.2 1. be fulfilled. For any e > 0 we can find 7 * and 
such that if 7 := J2x ^ox < 7 * there exists a spatial scale Lg and a temporal scale Lt, 
a set of good configurations Gp and a set of good Poisson process realisations Gu C 
25B4iAo)x{i,2,....Lt} P( (^Ibiock 4 (o 0 ) ^ ^u) > 1 — such that the contraction and 
coupling conditions ( |3.4| l, l |3.5| l from Sectionare fulfilled. Furthermore the random 
walk defined in l |4.11| l satisfies l |3.12| l in Section 


Proof. The crucial idea is that using the flow version ( |4.4| l we can augment the coupling 
argument in Lemma 13 in © to work with a set of (good) initial conditions 

:iGl} = {r]G : aj^ < r]{x) < /3/7 for x G B2l^(0)} (4.21) 

with a, (3 from ( |4.17| l and the (uncountable) index set I being defined implicitly here. 

The proof consists of 6 steps. For parameters 3> ATs > K't to be suitably tuned 
below, we set 

Ls = rA:slog(l/7)l 

Lt=L[ + L'l with L[ = \K[ log(l/7)l, L” = log(l/7)l • 


In the first step, we use the propagation properties of the deterministic system as de¬ 
scribed in Lemma [4^ together with the fact that for small 7 , the relative fluctuations of 
the driving Poisson processes are typically small to ensure that after time L[, the 'good 
region' has increased sufficiently. 

In the second step we use the flow version l |4.4| l and its contraction properties to 
ensure that in a subregion, after L" steps, coupling has occurred with high probability. 

Several copies of such subregions are then glued together in Steps 3 and 4. In 
Step 5 we use the fact that in a good region, the relative fluctuations of p are small so 
that Pp{k;x,y) is close to the deterministic kernel pxy; this ensures l |3.12| l. Finally, in 
the last step we collect the requirements on the various constants that occurred before 
and verify that they can be fulfilled consistently. 

Step 1. Let 


Xi := < max sup 

1 ||a;||.||y||<5LB,P^y>0,0<n<L; u>ao/~/ 


Ut^\u) 

Pxy'^ 


< 
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with tto from Lemma |4.4[ By standard large deviation estimates for Poisson processes, 
we have 


P(A’i) > 1 - {lOL^Rp^L^ exp(-cao/7) (4-22) 

(for some fixed constant c > 0 ) which can be made arbitrarily close to 1 by choosing 7 
small. 

By iterating l |4.17| l in combination with ( |4.18| l we see that 

A'l n ( 770 ( 0 ;) e [a/ 7 , 73 / 7 ] for a; G B 2 lS^)} C {r]L'^{y) G [a/ 7 ,/ 3 / 7 ] for j/ G S 5 l^( 0 )} 

(4.23) 

if the ratio L[/Ls is chosen sufficiently large. To verify this note that we can consider 77 
as a perturbation of the deterministic system C, from ( |4.14| l and on Xi the relative size 
of the perturbation is small when 7 is small (cf. @ Eg. (13) and the proof of Lemma 7]). 

Step 2. Let Go C I 3 ® 37 ,s(o)x{i,...,l"} j-j^g gg^. gj Poisson process path configurations in 
the space-time box f 33 L„( 0 ) x { 1 ,..., L/} with the property 


^ ( 4 >i,L"( 77 o))(a;) = ( 4 >i,L»(? 7 ''®^))(a;) for ||a;|| < Ls (4.24) 

with as in l |4.9| l and 77 ''®^ = 

Observe that for x G n G 


sup 77 ^ (x) - inf 77 W (a;) = sup ^ {f{y; - inf ^ {f{y; 

T • - T ' ^ iC. I ' ^ 


IGI 


iGl 


iGl - 


< ^ ( sup (/(y; 7 ?i*li)) - inf (/(y; 


iGl 


(4.25) 


^n-f ( sup fiy; ( jnf fir, Vn-i)) ■ 


Thus, 


E 


sup 77 , 
iei 


W(x) -inf77W(x) 

iGl 






sup /(y; 77 )/li) - inf f{y; 

lei 


(4.26) 


and we can now use contraction properties of / near its fixed point, analogous to the 
proof of (44), (45) in Lemma 13 in |j 6 |l. Note that if all r]‘^}_i{y), i G I are (locally around 
x) in the neighbourhood [a, /3] of m* (as reguired for ( |4.19| l), there are fjy^z G [a, / 3 ]^'*xG) 
such that 


sup f{y; 77 (;li) - inf f{y -^ 




I driiz] 


■/(y;yy .2 


(y) 



(4.27) 


Put 


V'rE(y) — l{r)(x)G[Q/7,/3/7] for ||a;||<B}- 

We have on {inf,g/ 7 /fl+H^+fl^(? 7 (;li) = 1} 

^ ® (sup7/i^(77W)77W(x)-inf7/fl(77W)77W(x) 


|Sb( 0 )| 


xeBii(o) 


n — l 


(4.28) 
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< 
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E[sup 77 ^ (a;) - inf ryW (a;) 


< 


l^fl( 0 )| i^i 

' ^ ’' xGBi?(0) ^ 

1^ Lni Y Y Py^ Y l^^/(y;’?^)l(sup??i-i(2)-inf77^1(0)) 

' x^BR(0)y^BR^{x) zeBa^iy) 

Y (sup7?ili(z) - inf 77ili(z)) Y Y P^v 


< 


2 G-Bh+Hp + h^(0 ) 

l^fl+-Rp+-RA( 0 )l 

|i?fi( 0 )| 


V^Br^ ( 2 ) 




(- 1 ) 


1 


2 G±Jfi+i?p + fi_ 5 ^ (U) 

XI (sup riYi.^) - ^i-i(^)). 


X (sup 77 i*li(z) - inf 77 ili(z)) 


(4.29) 


2 Gi^i?+Rp+K;^ (U) 

where we used < |4.26| l and ) |4.27| > in the second inequality and assume that R is so large 
that 

I(0) I C\ ^ \ ^ 1 /yi onx 

|B«( 0 )| - 2 ) ^ ^ *■ U.30) 

Note that the factor (1 — |) comes from l |4.19| l. 

We can iterate ( |4.29P for n = L", L" — 1,..., 1 to obtain on Xi (which in particular 
implies ^L,+k{Rp+R^){Vn-k) = 1 for /c = 1 , 2 ,..., tt - 1 ) that 


1 


|i?L.( 0 )| 


X E [ (sup 2jjL, (ryW, {x) - inf {rill )rl-ll {x)\ \Fq 

LXnfZT * * * * / J 

x&BlAO) 

u 


< c{e) 

< c(e) 


X (supr7«(z)-infr7«(z)) 

^iG/ * 6 / / 


7 

On the event {infig/ 4’r+l'^' (Rp+Rx)i'ilo'’) = 1} via Markov inequality l |4-31| l yields 
IP( max (supV'L,(77yO’7L((a;) - mfV'Ls(’7L()’7L((a;)) > 1 -^o) < - 

\||ai||<Ls i£j ^ ^ 4 ^/ t t / 

and the right hand side can be made as small as we like by choosing 7 small. 
Hence for dfa := { C/|b 3 ,.jo)x{i.....l"} ^ ^'o} we have 


(4.31) 


V{X2)>l-\BLMHef^' 


’ — a 
7 


) — a 
7 

(4.32) 


(4.33) 


Step 3. Let X 2 (i/, k) be the event that occurs in the space-time box whose 'bottom' is 
centred at (y,fc), i.e. X 2 {y,k) = {( 7 |B 3 p^^(j^)x{fe+i....,fc+L"} G Go}. By construction, on 




:=A’in f| X2{jL,ek,L[) 


(4.34) 


iG{-2.-l,...,2}. 

k—l,...,d 


we have 

r]o{x) G [ 0 / 7 ,/ 3 / 7 ] for X G H 2 l,( 0 ) ^ r]L,{y) = (4>i,l,( y) for || 7 /|| < 3Ls, 

(4.35) 
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i.e. I |3.4| ) holds. In particular 


P(-^3) > P(-^l) - 5'^(1 - P(A’2)) 


(4.36) 


Step 4. On 

Xi-.= X 3 n f| M0,jL''), (4.37) 

( |3.5| l holds as well, and we have 

P(A’4) > P(d:’3) - \L'JL”] (1 - P(A’2)) (4.38) 


Step 5. Note that on X 4 we have 

77o{x) G [a/7,^/7] for x G B2lA^) r]n{y) G [a/7, (3/7] for |j?/|| < 2L^,n= 1 ,... ,Lt. 

Then ( |4.10| l implies 

o. 13 

-j^Pxy <Pr,(k]x,y) < -pxy forx,y G B2l,( 0), fc = 


hence the total variation distance between x, •) and Px,. is at most (1 ~ V (£ — 1) 
uniformly inside this space-time block. We use Lemma [4.41 (in) to make this so small 
that coupling arguments as in the proof of Lemma |2.16| (with a comparison random walk 
that has a deterministic drift dmax ^ Ls/{L[ + L") show ( |3.12| l. 

Step 6. Finally, we verify that the constants Ks,KAK” can be chosen consistently so 
that all intermediate requirements are fulfilled. 


1. The right-hand side of l |4.22| l can be chosen arbitrarily close to 1 for any choice of 
Ks, K” by making 7 small. 

2. ( |4.23| l requires that K[ > j^Ks with sq from ( |4.17| l, ( |4.18| l. 

3. ( |4.31| l, which uses l |4.3Q| l L/ times, requires that Lg — (^Rp + R\)L” is large. This is 
achieved when Kg » {Rp + R\)K” (and 7 is small). 

4. The right-hand side of ( |4.33| l can be made close to 1 if K" > (—logc(e))“^ (and 
7 is small). This also implies that the right-hand side of ( |4.36| l can be chosen 
arbitrarily close to 1. 

5. For l |4.38| l note that ^ ^ is a fixed ratio when 7 is small, and 1 — P(A’ 2 ) can be 

made small by choosing 7 small. 


We see that for 7 < 7* for some 7* > 0, all requirements can be fulfilled e.g. by choosing 
iT/ := 2/(— logc(e)). Kg := C{Rp + R\)K” with some large C and K[ := j^Kg. □ 


5 Discussion of further classes of population models 

While we analysed in this article only one explicit spatial population model, namely lo¬ 
gistic branching random walks (LBRW) defined in l |4.4| l with the dynamics of the space- 
time embedding of an ancestral lineage given by ( |4.11| l, we do believe that the same 
program can be carried out for many related population models and that LBRW is in 
this sense prototypical; see also HE Remark 5]. We list and discuss some of these in the 
following five paragraphs. Note that implementing the details to verify the conditions 
of Theorem |3.1| for these models will still require quite some technical work and we 
defer this to future research. 
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More general 'regulation functions' The logistic function x i-)- x{m — Xx) (with A 
small enough) whose 'spatial version’ is used in the definition l |4.1| l can be replaced by 
some function (p : K+ —>• [0, a], a G (0, oo) U {00} with > 1, \\ mx ^ a ( t >{ x ) = 0 that 
possesses a unigue attracting fixed point x* = > 0. This will ensure that a result 

analogous to Lemma |4.4| holds, and an analogue of Proposition |4.7| can be obtained 
when a suitable small parameter is introduced. 

For example, in the ecology literature, in addition to the logistic model, also the 
Ricker model corresponding to cj){x) = xexp{r — Xx) and the Hassel model correspond¬ 
ing to ^{x) = mx/{l + Xx)^ are used to describe population dynamics under limited 
resources (r > 0, resp., m > 1 and 6 > 0 are parameters). Note that in all these cases, 
1 /A is related to a carrying capacity, so assuming A small means weak competition. 


More general families of offspring distributions As described in the informal dis¬ 
cussion above l |4.2| l, ( |4.5| l can be interpreted as stipulating that each individual at y in 
generation n has a Poisson number of offspring (with mean f {y; r]n)/Vniy)) which then 
independently take a random walk step. One could replace the Poisson distribution by 
some another family of distributions L{X(v)) on Ng that is parametrised by the mean 
E[2f(z/)] = v G [0,P] where D > supy f{y; y)/ri{y) and then define the model accord¬ 
ingly. If the family of offspring laws satisfies a suitably quantitative version of the law 
of large numbers (cf. Step 1 of the proof of Proposition |4.7p , one can derive an analogue 
of Proposition |4.7| 

For example, one could take an Ng-valued random variable X with mean E[X] = P 
and E[e“^] < 00 for some a > 0 and then define X{i') via independent thinning, i.e. 

= Y.f=i ^{Ui<u/v} where Ui, U 2 , ■ ■ ■ are i.i.d. uniform([0,1]). 


'Moderately' small competition parameters As it stands. Theorem |4.3| requires 
sufficiently (in fact, very) small competition parameters (cf. Assumption |4.2[ l. This is 
owed to the fact that our abstract 'work-horse' Theorem |3.1 [ requires (very) small eu in 
Assumption |3.3| and e in Assumption |3.9| (we in fact did not spell out explicit bounds). 
In simulations of LBRW one observes also for moderately small competition parameters 


Xxy apparent stabilisation to a non-trivial 'equilibrium' as required by Assumption 3.3 


We note that the assumptions of Theorem [37^ are 'effective' in the sense that they only 
require controlling the system (rjn) and the walk in certain finite space-time boxes. 
Thus, a suitably quantified version of Theorem |3.1| allows at least in principle to as¬ 
certain via simulations that for a given choice of parameters m, (pxy) and {X^y) the 
system (? 7 „) has a unique non-trivial ergodic equilibrium and that the conclusions of 
Theorem 14.31 hold. 


Continuous-time and continuous-mass models An infinite system of interacting 
diffusions that can be obtained as a time- and mass-rescaling of LBRW is considered in 
nil, see Definition 1.3 there; one can in principle define an 'ancestral lineage' in such a 
model which will be a certain continuous-time random walk in (the time-reversal of this) 
random environment. It is conceivable that a coarse-graining construction similar to the 
one discussed here can be implemented and that in fact an analogue of Theorem |4.3| 
can be proved at least for suitably small interaction parameters. 

Reversible Markov systems Obviously, we tailored Theorem [3.1 [ and its assumptions 
to a random walk that moves in the time-reversal of a non-reversible Markov system rj 
which possesses two distinct ergodic equilibria, our prime example being the (discrete 
time) contact process. 
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If we instead assume that ry is a reversible Markov system with local dynamics and 
'good' mixing properties possessing a unique equilibrium (for example, the stochastic 
Ising model at high temperature), the assumptions from Section |3.1| will be fulfilled 
as well - this case is in fact easier since 'good blocks' in 77 (k+i)Lj as required in ( |3.4| l 
of Assumption |3.3| will have uniformly high probability anyway, irrespective of rjuLi ■ 


Assume in addition that we can verify Assumption |3.9| for the walk. For example, this 
can be done by requiring that the walk is a (sufficiently small) perturbation of a fixed 
symmetric random walk or by assuming a (small) a priori bound on the drift. Then we 
don't need to require the symmetry assumption ( |3.13| l (in fact, the resulting walk can 
have non-zero speed). This re-reading of Theorem |3.1| and its proof allows to recover a 
special case of 11271 Thm. 3.6] where, using entirely different methods, a CLT is obtained 
for random walks in dynamic environments that satisfy sufficiently strong coupling and 
mixing properties. 


A An auxiliary result 

The following result should be standard, we give here a brief argument for complete¬ 
ness' sake and for lack of a precise point reference. 

Lemma A.l. Let T = (7^n)n=o.i,... be a filtration, T, T' finite J^-stopping times, and Y a 
bounded random variable. We have 


E [Y I J~t\ ll{r=T'} = E [Y \ J~t'\ ll{r=T'} a.s. (A.l) 

Proof. Note that {T = T'} e J't H Tt’ because 
{T = T'} n{T = n} = {T = T'} n {T' = n} = {T = n} n {T' = n} G Tn, n = 0,1,.... 
Furthermore we have 

A € j-T u J-T' A n {r = T^} G j-T n j-t' . 

To see this note that for A e Tt 

An{T = T'}n {T' = n} = {An{T = n}) n {T' = n} G n = 0,1, ■ • ■ ■ 

Thus, we obtain A n {T = T'} G and a similar argument for the other case shows 
the assertion. By approximation arguments we find that 

Z is J'T-measurable Z1^x=t'} is (^t H Jr'l-measurable. 

Let Z be a version of E I ^t\ = 1 {t=t'}E[F | Ft], i-e., Z is Jr-measurable, 

E [Z1a\ = E for all A G Ft- We may assume that Z = Z1^t=t'}- Then, Z 

is also aversion of E | Ft'] = ] 1 {t=t'}E[F | Ft']- Furthermore Z = Z 1 {t=t’} 

is also Jt’'- measurable and for A' G Ft', 

E [Zl,4/] = E \ZlA'n{T=T'}\ E [i'^l{T=T'}lA'n{T=T'}] = E ■ 

This concludes the proof of the lemma. □ 
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